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Abstract. One of the principal topics of this paper concerns the realization of 
self-adjoint operators Lq q in L-^ [Q; d" x)"^ , m, n € N, associated with diver- 
gence form elliptic partial differential expressions L with (nonlocal) Robin-type 
boundary conditions in bounded Lipschitz domains Q C M". In particular, we 
develop the theory in the vector- valued case and hence focus on matrix- valued 
differential expressions L which act as 

= ~{ ^ ^A^^'tk^kUp]] , u = {ui,. . . ,U,n). 

The (nonlocal) Robin-type boundary conditions arc then of the form 

[{du/du) + eu]\Qn=0, 

where Q represents an appropriate operator acting on Sobolev spaces asso- 
ciated with the boundary dQ of Q, and d/df denotes the outward pointing 
normal derivative on dil. 

Assuming ^ in the scalar case m = 1, we prove Gaussian heat kernel 
bounds for Lq q by employing positivity preserving arguments for the asso- 
ciated semigroups and reducing the problem to the corresponding Gaussian 
heat kernel bounds for the case of Neumann boundary conditions on dQ. We 
also discuss additional zero-order potential coefficients V and hence operators 
corresponding to the form sum Lq q + V. 
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1. Introduction 



In a nutshell, this paper is devoted to a new class of self-adjoint realizations 
in (Px) of elliptic partial differential expressions in divergence form, 



on bounded Lipschitz domains C R", n ^ 2, with Robin boundary conditions 
of the form [{du/di^) + Ou]\dii = 0. (Here denotes the outward pointing imit 
vector and ^ is a suitable function on the boundary dfl of Cl.) Following [69], 
we put particular emphasis on developing a theory of nonlocal Robin boundary 
conditions where the function 9 on dfl is replaced by a suitable operator Q acting 
in L^{dfl;d"'~^oj), with d"'~^u) representing the surface measure on dfl. (More 
precisely, Q acts in appropriate Sobolev spaces on the boundary of 17, cf. Section 
3.) The resulting self-adjoint operator in L'^{Q;(Px) is then denoted by io,a and 
we study its resolvent and semigroup, proving a Gaussian heat kernel bound and 
a corresponding bound for the Green's function of Lq q on the basis of positivity 
preserving arguments. 

The corresponding case of Dirichlet boundary conditions, and similarly (al- 
though, to a somewhat lesser extent), the one of Neumann boundary conditions 
on 9f2 has been extensively studied in the literature. An authoritative survey of 
this area until about 1989 was written by Davies [50], and newer developments 
since then were treated by Ouhabaz [119] in 2005. For the study of elliptic opera- 
tors on Lie groups and the associated semigroup kernels we refer to the monograph 
by Robinson [127]. A thorough study of the heat equation and the heat kernel of 
the Laplacian on Ricmannian manifolds has been undertaken by Grigor'yan [79], 
and very recently, Dirichlet and Neumann heat kernels and associated two-sided 
bounds were developed for inner uniform domains by Gyrya and Saloff-Coste [85]. 

The case of Robin-type boundary conditions on the other hand, received much 
less attention. The latter experienced a boost due to the seminal paper by Arendt 
and ter Elst [15] in 1997, and the past fifteen years have seen a number of interesting 
developments in this area. 

Before specializing to matters related to heat kernel bounds, we digress a bit and 
mention some of the literature devoted to Robin boundary conditions: A probabilis- 
tic approach to Robin- type (or third) boundary value problems for L = A + fe-V + g 
was undertaken by Papanicolaou [121] for bounded C'^-domains ft; the case of the 
Robin problem for the Laplacian on general domains, including fractals, was treated 
by Bass, Burdzy, and Chen [25]. The short-time asymptotics of the heat kernel for 
the Laplacian with Robin boundary conditions on a smooth domain was discussed 
by Zayed [148] . Inequalities between Robin and Dirichlet eigenvalues for the Lapla- 
cian on bounded Lipschitz domains Q were derived by Filonov [61]. (The paper by 
Filonov motivated two of us to extend such results to the case of nonlocal Robin 
Laplacians which lead to their introduction in [69].) An isoperimetric inequality 
similar to the Fabcr- Krahn inequality for Robin Laplacians (and more generally, in 
the case of mixed Dirichlet and Robin boundary conditions on parts of the bound- 
ary) on bounded Lipschitz domains was established by Daners [41]; in this isoperi- 
metric context we also mention work by Kennedy [96], [97], [98], [99]. Uniqueness 
in the Faber-Krahn inequality for Robin Laplacians was derived by Daners and 
Kennedy [43] , and for a recent alternative approach to the Faber-Krahn inequality 
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for p-Laplacians with Robin boundary conditions we refer to Bucur and Daners 
[30]. Domain perturbations for elliptic equations with Robin boundary conditions 
were studied by Dancer and Daners [38]; domain monotonicity for the principal 
eigenvalue of Robin Laplacians was investigated by Giorgi and Smits [72], [73]. 
Spectral stability of Robin Laplacians and Laplacians with mixed boundary condi- 
tions on bounded domains (Lipschitz and more general than that) defined in terms 
of quadratic forms were treated by Barbatis, Burenkov, Lamberti, and Lanza de 
Cristoforis in various collaborations [24], [32], [33]. Eigenvalue asymptotics for 
Robin Laplacians was studied by Daners and Kennedy [44]. Schatten-von Neu- 
mann properties of resolvent differences and the asymptotic behavior of singular 
values of nonlocal Robin Laplacians on bounded smooth domains were investigated 
by Behrndt, Langer, Lobanov, Lotoreichik, and Popov [26]; in the case of local 
Robin boundary conditions these results were subsequently improved and extended 
to strongly elliptic symmetric partial differential operators by Grubb [82] (see also 
[83]). Nodal line domains for Laplacians with (local) Robin boundary conditions 
on bounded Lipschitz domains were studied in [100]. The Robin boundary value 
problem for the two-dimensional Laplacian in domains with cusps was studied by 
Kamotski and Maz'ya [92]. Hardy inequalities for Robin Laplacians were derived 
by Kovafik and Laptev [104]. The problem of minimizing the nth eigenvalue of the 
Robin Laplacian was recently investigated by Antunes, Freitas, and Kennedy [7]. 

Returning to Gaussian heat kernel bounds for uniformly elliptic partial difi^er- 
ential operators, the literature associated with Dirichlet boundary conditions is 
so enormous that one would be hard-pressed to do it justice in the introduction 
to a manuscript of the present scale. Instead, we refer to [50] and [119] which 
nicely survey the state of literature until about 1989 and 2005, respectively. Thus, 
we will now focus on studies in connection with Robin boundary conditions, the 
principal topic of this paper. We already mentioned the influential paper by 
Arendt and ter Elst [15] on Gaussian heat kernel bounds for L^.n on bounded 
Lipschitz domains with local Robin boundary conditions indexed by the function 
9 on d^l. The authors study the case of real, not necessarily symmetric matrices 
aj,k G L°°(ri; d^x), 1 ^ j,k ^ n, satisfying a uniform ellipticity condition, assuming 
^ 9 £ L°°{dQ;d"'~^uj), and permitting lower-order (complex-valued) coefficients 
under some smoothness hypotheses. Subsequently, Daners [40], in the case of real- 
valued coeflicients, extended their results by removing the smoothness assumptions 
on the lower-order coefficients. Moreover, Ouhabaz [118] and [119, Sect. 4.1, Ch. 6], 
further extended these results by permitting complex-valued coeflicients aj^k in L 
which introduces a variety of challenges. One should also mention that the results 
developed by these authors also permit mixed boundary conditions in the sense 
that one may have Dirichlet boundary conditions at a part of the boundary and 
Neumann or Robin boundary conditions on the rest of the boundary. In addition, 
local Robin-type problems were studied by Daners [39] on arbitrary bounded and 
open sets Q C M"; he proved positivity improving of the underlying semigroups if 
fl is connected and < ^ e L°°(9f2; d"~^a;). The case of arbitrary bounded and 
open sets H, C M" was also treated by Arendt and Warma [16]. [17], [141] (sec also 
[142]), who developed a general theory with 9 replaced by a (positive) measure on 
the boundary, and who proved Gaussian heat kernel bounds in this context. In par- 
ticular, a characterization of generalized local Robin Laplacians whose semigroups 
satisfy domination by the Neumann semigroup (cf. (4.21) for details) was provided 
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in [17]. Positivity preserving semigroups in the case of bounded Lipschitz domains 
fl and sign indefinite 6 G L°°{dCl;(P~'^oj) were established by Daners [42] (and 
again mixed boundary conditions were permitted). For an illuminating survey on 
semigroups and heat kernel estimates we also refer to [11]. 

Next, we briefly recall some facts concerning positivity preserving integral opera- 
tors. For simplicity, we describe the concrete case of linear operators in the Hilbert 
space 1/^(0; rf" a;) (f2 c M" a bounded Lipschitz domain) and refer to Section 2 
for the abstract setting: A bounded operator A in L^{fl;d"'x) is called positivity 
preserving if 

O^fe L^{n; d"x)\{0} implies Af ^ a.e. (1.2) 

This will be denoted by A :^ (resp., by =^ A). More generally, A )^ B (resp., 
B =^ A) then abbreviates that ^ — B is positivity preserving. A key result in this 
context then concerns the fact that a bounded integral operator A in L^(J7;d"a;) 
is positivity preserving if and only if its integral kernel A{-, •) is nonnegative a.e. 
on f2 X O. In particular, if A and B are bounded integral operators in L'^{Q.; dJ^x), 
then 

Q^B ^ A if and only if s$ B{-, ■) A{-, ■) a.e. on 17 x f7. (1.3) 

In particular, =^ i? =^ A yields a bound on the integral kernel of B in terms of 
that of A. 

Our approach to proving a Gaussian heat kernel bound for the class of operators 
I/e,o then rests on the following strategy: Assuming 6 ^ 0, we will prove that 
the semigroups e"*^'^ ", t ^ 0, are positivity preserving, and so are the differences 
g-tie.n _ g-tiN,s7^ t ^ 0, where L^^a denotes the case of Neumann boundary 
conditions (i.e., the special case 9 = 0). Thus, assuming ^ 0, we will prove 

=^ e-*-^«-" e"'-^-^'-", t > 0, (1.4) 

and hence obtain the Gaussian heat kernel bound for Le.n by means of (1.3) and 
the corresponding known Gaussian heat kernel bounds for Ln^q. 

We briefly turn to a description of the contents of each section: Section 2 is 
devoted to an abstract discussion of positivity preserving (and improving) integral 
operators with special emphasis on resolvents and semigroups. Much of the material 
in this section revolves about an important early paper by Davies [46] and the 
remarkable papers by Bratteli, Kishimoto, and Robinson [28] and Kishimoto and 
Robinson [101]. 

Section 3 develops in detail the theory of nonlocal Robin-type operators Le.n 
on bounded Lipschitz domains c R". In fact, we will go a step further and treat 
the m X m matrix-valued case, m e N, where L acts like 

Lu = -( ^ ( E ^j,k^kU0 ]] , u = {ui,...,Um)- (1.5) 

While we will naturally follow the outline of an earlier treatment in the special scalar 
case m = 1 of the Laplacian, L = — A in [69], we emphasize that the presence of 
the tensor coefficient A = {a"f) i^j,k^n € L°°{fl;d"'x)"^^'^ requires a careful re- 

examination and extension of the earlier arguments in [69] . We also note that our 
results are of interest in the special case of local Robin boundary conditions where 
9 corresponds to the operator of multiplication by the function 6 defined on 
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dft as we do not have to assume that 6 G L°°{d^;d"^^u!), but arc able to permit 
appropriate LP-conditions on 0. 

In our final Section 4 we then derive the Gaussian heat kernel bounds for Le.o 
employing the positivity preserving (in fact, positivity improving) arguments as 
outlined in connection with (1.4). This section concludes with a scries of remarks 
that enhances the main result on Gaussian heat kernel bounds, including the ad- 
dition of a nonnegative potential term ^ V G Ll^^{Q; (T" x) to Lq^q using the 
method of forms, and appropriate negative and small form potential perturbations. 

Appendix A summarizes the principal results of Sobolev spaces on Lipschitz 
domains O c and on their boundaries dfl, and Appendix B is devoted to 
the basics of scsquilinear forms and their associated operators. Both appendices 
provide the basics for Section 3 and are included to guarantee a certain degree of 
self-containment of this paper. Appendix C recalls some bounds for heat kernels 
and Green's functions; in particular, we provide a streamlined approach to Green's 
function bounds including the case of dimension n = 2. 

Finally, we briefly summarize some of the notation used in this paper: Let 
be a separable complex Hilbert space, (•, •)-« the scalar product in H (linear in the 
second argument), and I-h the identity operator in H. 

Next, if T is a linear operator mapping (a subspace of) a Hilbert space into 
another, then dom(T) and ker(T) denote the domain and kernel (i.e., null space) of 
T. The closure of a closable operator S is denoted by S. The spectrum, essential 
spectrum, discrete spectrum, and resolvent set of a closed linear operator in a 
Hilbert space will be denoted by a{-), fTcss(-): '^di'), and p(-), respectively. 

The convergence in the strong operator topology (i.e., pointwisc limits) will be 
denoted by s-lim. Similarly, limits in the weak (resp., norm) topology are abbrevi- 
ated by w-lim (resp., by n-lim). 

The Banach spaces of bounded and compact linear operators on a separable 
complex Hilbert space are denoted by B{'H) and i3oo('H), respectively; the cor- 
responding ^^'-based trace ideals will be denoted by Bp{'H), p > 0. The trace of 
trace class operators in H is denoted by tifi{-). The analogous notation 6(^*1, A2), 
Bao{Xi,X2), etc., will be used for bounded and compact operators between two 
Banach spaces Xi and X2. Moreover, Xi ^ X2 denotes the continuous embedding 
of the Banach space Xi into the Banach space A2. 

Given a u-finite measure space, (M, A^, /x), the product measure on M x M will 
be denoted by /x ® /U. Without loss of generality, we also denote the completion 
of the product measure space (M x M,M. ® M, ® fj) by the same symbol and 
always work with this completion in the following. For brevity, the identity operator 
in i^(Af;rf/i) is denoted by /a/, the scalar product and norm in L?{M]dii) are 
abbreviated by ( • , • )m and || • ||m. whenever the underlying measure is understood. 
For a linear subspace T) of L'^ {AI ; dn) , the cone of nonnegative elements in V is 
denoted by D+. Lastly, Xs denotes the characteristic function of the set S. 

2. On Positivity Preserving/Improving Integral Operators 

We start by recalling some basic facts on positivity preserving/ improving oper- 
ators. Throughout this paper we will make the following assumption. 

Hypothesis 2.1. Let {M,M.,fi) denote a a-finite, separable measure space asso- 
ciated with a nontrivial measure {i.e., < /x(M) < 00). 
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Assuming Hypothesis 2.1, L'^{M; djj) represents the associated complex, separa- 
ble Hilbert space (cf. [27, Sect. 1.5] and [91, p. 262-263] for additional facts in this 
context). 

The set of nonnegative elements ^ / e L'^{M;diJ,) (i.e., f{x) ^ /U-a.e.) is a 
cone in L^{M; dfi), closed in the norm and weak topologies. 

Definition 2.2. Let A he a linear operator A in L'^{M;dfi). Then A is called 

positivity preserving {resp., positivity improving) if 

7^ / e dom(A), / ^ n-a.e. implies Af ^ {resp., Af > 0) n-a.e. (2.1) 

The operator A is called indecomposable if {0} and L^{M;dfj,) are the only closed 
subspaces of L'^{M;d^) left invariant by A. 

Positivity preserving (resp., improving) of A will be denoted by 

A > (resp., A >- 0). (2.2) 

(or by =^ ^ (resp., -< A)). More generally, ii A,B G B{L'^{M; d/i)), then 

A^ B (resp.. Ay B) (2.3) 

(or B ^ A (resp., B -< A)), by definition, imply that A — B is positivity preserving 
(resp., positivity improving). We note parenthetically, that in large parts of the 
pertinent literature, "positivity preserving operators" are just called "positive", 
but due to the obvious conflict of notation with the concept of positive operators in 
a spectral theoretic context, we follow the mathematical physics community which 
prefers the notion of positivity preserving. 

Moreover, if A has the property that whenever 7^ / G dom(j4) then also 
I/I G dom(A) (in particular, if dom(j4) — L'^{M; dji)) one infers that 

A ^ if and only if \Af \ ^ A\f \ /x-a.e. (2.4) 

(upon decomposing / = /+ — /_, f± = (|/|±/)/2). Similarly, one concludes that 

^ ^ (resp., A >- 0) if and only if {f,Ag)M > (resp., {f,Ag)M > 0) 

for all < / e i^(M; d/i)\{0}, ^ 5 e dom(^)\{0}, ^'^'^^ 

using the fact that 

g e Z/^(M; rf/u)\{0} is nonnegative (resp., strictly positive) if and only if 

(/, 9)m > (resp., (/, g)M > 0) for all < / e ^^(M; dfx)\{Q]. ^^'^^ 

This, in turn, follows from cr-finitcness of M, that is, M = UneN^^"' where Af„ G 
M, n{Mn) < 00, defining S = {x G M \ g{x) <Q), Sn = S Mn, n G N, and from 
the fact that if n{Sno) for some no G N (recalling that /x(M) > 0), 



( resp., < ) (xs„„ , 5)m = / Xs„„ (2^)5(2;) dn{x) 0. (2.7) 

JM 

Turning our attention to integral operators in LF'{M; dfi) with associated integral 
kernels A(-, •) on M x M, we assume that 

A{-, •):MxM^Cis/i0 ;U-measurable, (2.8) 

and introduce the integral operator A generated by the integral kernel A{-,-) as 
follows: 

{Af){x) = I A{x,y)f{y)dn{y) for /x-a.e. x€M, 
Jm 
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/ e dom(A) = \gGL^{M; dji) 



A{x, ■)g{-) is integrable over M for /^-a.e. x G M, 

and/ A{;y)g{y)dfx{y)eL\M;dA. (2.9) 
Jm ) 

As shown in [91, Theorem 11.1], 

AeB{L^{M;d^i)) whenever doni(A) ^^(M; d^) (2.10) 

(as A turns out to be closed in this case). In addition, we also note that by [91, 
Theorem 11.2], |A(-,-)] generates a bounded integral operator in L'^{M;dn) if and 
only if 

/ \f{x)A{x,y)g{y)\d{ti ® ti){x,y) <oo, f,g G L\M;dn). (2.11) 

J MxM 

For additional resuhs on integral operators we refer, for instance, to [56, Sect. 
9.5], [86], [103], [143, Ch. 6]. 

Next, we state the foUowing resuh (special cases of which appear to be well- 
known, but we know of no published proof at this instant). 

Theorem 2.3. Assume Hypothesis 2.1 and suppose {without loss of generality) 
that M can be written as M = U^gj^M„, where Mn G A4, ^{Mn) < oo, n G N. 
Moreover, suppose that A is an integral operator in L'^{M;dfj,) generated by the 
integral kernel A{-,-) on M x M satisfying 

Ai-,-)\M^xM„ G L\MmX Mn,Mm<S^Mn;d{lJ,ra<S) IJ-n)), m,n GN, (2.12) 

where 

Mn = M„nM = {Mnr]S\SGM} = {TCMn\TGM}, n e N, (2.13) 

Atn=M|M„, nGN. (2.14) 
Then A is positivity preserving if and only if 

A{-,-)^ IX® ix-a.e. on M X M. (2.15) 

Proof. Clearly, a positivity preserving integral operator of the type (2.9) must have 
a real- valued integral kernel A{-, ■). SuiEciency of the condition (2.15) for positivity 
preserving of A then follows directly from the representation of A in (2.9). 

In order to prove necessity of the condition (2.15) it suffices to proof that posi- 
tivity preserving of A implies that 

^(•,-)|m„xm„ > (g)/i„-a.e. on M„ X M„, m,neN. (2-16) 

Next, for each (m, n) G N^, we introduce the set 

Sm,n = {{x, y)GMmXMn\ A{x, t/) < Mm O /U„-a.e.}. (2.17) 

Then Sm,n G Mm ® Mn, (m, n) G N^. 

By [91, Lemma 11.1.2], for each e > 0, there exist disjoint measurable rectangles 
Sj{s) X Tj{e) G Mm-^Mn, with iimiSj{s)) < oo, UniTjie)) < oo, j = 1, - •• ,N{e), 
N{e) G N, such that with 

N(s) 

Rme) = U [^^■(^) ^:'(^)]' (2-18) 
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one infers that 

/ A{x,y)d{iim® li'n){x,y) - I A{x,y) d{nm Hn){x,y) ^ £• (2.19) 
Next we claim that 

/ A{x,y)d{iJm®fin){x,y)>0. (2.20) 

Indeed, since ^(•,-)|m™xm„ G L^{M,n^ Mn,Mm<^Mn;d{fim(^ fin)), a^^d Xs^(._) ^ 
L2(M„; dfim), Xtj^s) ^ L^iMn.dfin) as ,u„(S'-,(e)) < oo, iJn{Tj{e)) < oo by hypoth- 



esis, an application of Fubini's theorem yields 



/ A{x,y)d{iJ,m®IJ'n)ix,y) (2.21) 

JSiie)xTj{e) 

= {Xs,,,_,,Ax^^J^^O, J = !,■■■ ,Nie), (2.22) 

as A is assumed to be positivity preserving. Thus, combining (2.19) and (2.22) one 
concludes that 

/ A{x,y)d{iJ,m'® Hn){x,y) 

JSm.n. 

= (/ A{x,y)d{iim'® lin){x,y) - j A{x,y) d{ij,m 'Si iJ.n){x,y)] 
+ A{x,y)d{iJ.mS iJ.n){x,y) e[-e,oo). (2.23) 

Letting £ 4- 0, (2.23) finally yields that 

/ A{x,y)d{iJ.mSiJ.n){x,y)^0 (2.24) 

JSm.n 

(AtmOMn)(5'm,„) = 0. (2.25) 

(m8)m)( U Sm,nj=0, (2.26) 

implying (2.15). □ 

Remark 2.4. In connection with the sets M„, n G N, in Theorem 2.3, which arc 
used to formulate a substitute for the lack of local integrability of the integral kernel 
A{-,-) (due to the absence of any topology imposed on M), we note that one can 
assume that the Mn are mutually disjoint. Mm O Mn = 0, m 7^ n, or else, that they 
are nesting, M„ C M„_|_i, m,n € N. In addition, we note that upon introducing 
Ln = Uj=i -^j' n gN, the Cartesian product of M with itself takes on the simple 
form MxM = U„eN [^n x i„] . 



and hence 
Consequently, 
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Remark 2.5. It is clear that if the integral kernel for A satisfies A{-, •) > /i(8)/x-a.e., 
then A is indecomposable (cf. [47, Sect. 7.1]), in fact, positivity improving, that is, 
^ :^ 0. The converse, however, is clearly false as the following elementary example 
shows. 

Example 2.6. In the Hilbert space i^([0, l];dfio) with 

Mo(S)= / dx + XB(l/2), Be!B(R), (2.27) 
Jb 

where f8(K) denotes the Borel a-algebra on K., and Xsi') characteristic func- 
tion of the set S, consider the rank-one operator Aq with associated integral kernel 
Ao(-, •) given by 

Ao{x,y) = ao{x)ao{y), ao(t) = |t - (1/2)|, x,y,tG [0,1]. (2.28) 
Then Ao{-, •) ^ and Aq >- 0, but clearly Ao{-, •) is not strictly positive yUo^/io-a.e. 

Indeed, using the fact that h G L^{M ; d/x) is strictly positive (i.e., h > /z-a.e.) 
if and only if {f,h)M > for all < / e {M ; diJ,)\{0} , one concludes, upon 
identifying M = [0, 1] and dfx with Lebesgue measure dt, that 

/ ao{t)f{t)d^io{t)^ I ao{t)f{t)dt>0, < / e i'([0, 1]; dMo)\{0}, 

J[0,1] J[0,1] 

(2.29) 

implying 

(/,^o5)[o,i] = / ao{x)f{x)dfio{x) ao{y)g{y) di^o{y) 

J[0,1] J[0,1] 



aoix)f{x)dx aoiy)giy)dy>0, (2.30) 
[0,1] J[o,i] 

0</,ff eL2([o,l];d^o)\{0}. 

One observes that the subset of positivity preserving operators of B{7i), 

{AeB{n)\A)pO} (2.31) 

is a cone, closed under multiplication, taking adjoints, and under taking weak 
operator limits. 

As an immediate consequence of (2.5), we note that if A,B G B{H) then 
A ^ if and only if A* ^ 0, 

A ^ B ^ if and only if A* > B* > 0, (2.32) 
if A ^ 0, S > then AB > 0. 

Wc also recall the following (slight refinement of a) useful result discussed in [46] 
on domination in connection with trace class and Hilbert- Schmidt operators and 
we add one more item with respect to compactness : 

Lemma 2.7. Assume Hypothesis 2.1 and suppose that 

A,B& B{L'^{M; dfj,)) and A>p B >p= Q. (2.33) 

(i) The following norm bound holds, 

||-4||B(i,2(M;d^)) > \\B\\B(L-^(M-4ti.))- (2.34) 
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(ii) Assume in addition to (2.33) that A S B2n{L^{M; dn)) for some n G N. Then 
also B e B2n{L'^{M;dn)) and 

\\MB2n{L^{M;dij.)) > \\B\\B2n{L^{M;diJ.))- (2.35) 

(Hi) Assume in addition to (2.33) that A G Boo{L'^{M;dfi)) . Then also B e 

B^{L\M■d^i)). 

(iv) Assume in addition to (2.33) that A ^ 0, B ^ 0, and A e Si (L^(M; d/z)) . 
Then also B e Bi{L'^{M\d^i)) and 

tl"l,2(M;d/i)(^) = imiBi(Z,2(M;d/i)) ^ l|-B||Bi(Z,2(M;d^)) = L'^ (M;dn){B) ■ (2.36) 

Proof, (i) We refer to [46, Lemma 1.2] for the proof of (2.34) (see also the paragraph 
following Theorem 2.10 for a simple argument). 

(a) Repeatedly employing (2.32) one infers that A)^ B )pQ implies A* A )p B*B ^ 
and hence 

(A* A)" ^ (B*S)" :>= 0. (2.37) 
Recalling that T e -B2n('H) if and only if T*T G Bn{H), which in turn is equivalent 
to (T*T)" eBi(^) and that 

Wnt^iu) = \\T*m^^^^ = ||(T*T)"||g^(„) = tvn{{T*Tr), (2.38) 

one can apply item (iv) (with the pair A,B replaced by {A*A)'^, {B*BY) to con- 
clude that B G B2n{L\M; dy)) and 

t^L^iM-4,){{A*Ar) = ||(AM)"||g^(^,(^^^^jj 

Item (zm) is a special case of a result proven in [53] and [125] (see also [107]). 
(if) Following [46, Lemma 1.2], one constructs a filtering increasing set of orthog- 
onal projections Pl in L'^{M; dfj) associated with closed subspaces L generated by 
finitely many characteristic functions such that P^ strongly converges to Im ■ Since 
the details of this construction arc a bit involved, we now pause and present the 
precise argument: First, one recalls from Hypothesis 2.1 that {M,M.,^) is a a- 
finite measure space. Consider = {E G A1 | < /i(-E) < oo} and, for every set 
E G M*, define xe = I^{E)^^^^Xe' that is, the L-^-normalization of the character- 
istic function of E. Next, for every finite family {Ei, . . . ,E]y} of mutually disjoint 
sets in A^* introduce 

L{Ei,...,En) = (Y.''i^^^ ajGM, Ki^Tvj, (2.40) 

that is, the linear span of the xb/s. Going further, introduce 

Z{M,M,^i) = {L{Ei,...,En)\N e Nand{£;i, . . . C M,}, (2.41) 

and equip this family of sets with the partial order induced by the inclusion. Finally, 
for each L = L{E\, . . . , En) G £(M, A^, fj) denote by Pl the operator of projection 
in L'^iM-.d^j.) onto the linear (closed) subspace L of L^{M;diJ,). Hence, given that 
the collection {xB^ji^j^Ar is an orthonormal basis for L{Ei, . . . , En), for every 
function / G L^{M; dfj), one has 
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Moreover, for every Li,L2 G £,{M,M,ii) such that L2 C Li, the Pythagorean 
theorem implies 

11/ - PlJW'm = 11/ - PlJW'm + \\PlJ - PlJW'm (2.43) 
for every / e L'^{M; d/z). This imphes that for every / e L^{M; d^) the bound 

\\f-PLj\\M^\\f-PLj\\M, (2.44) 

holds whenever Li,L2 S £(M, A^,/z) are such that L2 C Li. At this stage, we 
make the claim that the filtering increasing family of projections {Pl} Les,{M,M,iJ,) 
converges strongly in L'^{M;dfj,) to the identity on M, that is, 

s-limPi=7M in L'^{M;dn). (2.45) 

To justify the above claim, fix an arbitrary function / S L'^{M; d^), along with 
an arbitrary number e > 0. Also, for each A e (0,oo) set fx :— Jxb^ where 
Ex ■■= {x & M \ \ f(x)\ ^ A}. Hence, Ex & M for every A € (0,oo). Since the 
sequence {fx} x>o converges to / in L?{M; d^) as A — >■ 00 (by Lebesgue's Dominated 
Convergence Theorem), one can find Aq € (0, 00) with the property that 

||/-/aoI|m<£/8. (2.46) 

Next, given that the measure space (M, A^,/z) is cr-finitc, there exists a family 

{AffejfeeN of mutually disjoint sets in A^, with the property that M — yj Mk- 

fceN 

Since the sequence {fxoXu" m }neN converges to fxo in L^{M;dii) as n — >■ 00 

fe = l 

(again, by Lebesgue's Dominated Convergence Theorem), one can find no G N such 
that 

ll/Ao-/AoX,no „ I|m<£/8. (2.47) 

To proceed, abbreviate Mq = X no and /o = fxoXu ■ Also, pick an arbitrary 

integer n such that n > Aq (whose actual value is to be specified later). In particular, 
|/o(a;)| < n for every x G M, and fo{x) = for every x G M\Mq. Also, from (2.46) 
and (2.47) one infers that 

||/-/o||m <e/4. (2.48) 
For each integer fc G [1 — n2",n2"], we now define 

i^„,fc = /o"'([^,^))nMo. (2.49) 
By design, for every A; G [1 — n2", n2"] fl Z we then obtain 

En,k e M, En,k C Mo, (2.50) 

and ^ < fo{x) < ^ for each x G £„,fe. (2.51) 

Furthermore, 

the family {En,k}i-n2'^^k^n2^ is a disjoint partition of Mq. (2.52) 

Hence, whenever k G [1 — 712", n2"]nZ is such that iJL{En,k) > 0, we have En^k € -M.* 
and 

^ < KEn,k)-' [ fodfx^ ^. (2.53) 
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In light of (2.51), this shows that whenever fc G [1 — n2",n2"] n Z is such that 
^{En,k) > we have 



^2 "for each x € En,k- 



(2.54) 



Next, thin out the family {iSn,fe}i-n2"^fe^n2" by throwing away the measure zero 
sets, and relabel the remaining ones as {Ex, . . . ,En}. In addition, consider L = 
L{Ei, . . .,En) e &{M,M,fi). Because of (2.52) it follows that {Ei,..., En} is, up 
to a set of measure zero, a disjoint partition of Mq. Based on this, the fact that /o 
vanishes outside Mq (as pointed out earlier), as well as (2.54) and (2.42), we may 
then estimate 

N 



\Mx) - {PLfo){x)\ 



E/o(a;)x./a:)- (^i^,)-! / fod^i 

j = l ^ •'^3 

N 



(2.55) 



for ii-a.c. X e M. Since both /o and Pl/o vanish identically outside Mq, we deduce 
from (2.55) that 

ll/o - Pl/oIIm < 2-XMo)^/^ (2.56) 

Consequently, 

11/ - PLfU < 11/ - Mm + ll/o - PlMW 
+ ||Pl(/o-/)||m 

<(e/2) + 2-XMo)i/2, (2.57) 

by (2.48) (used twice) and (2.56). Hence, if the integer n G (Ao,oo) is chosen 
large enough so that 2~"/x(Mo)^/^ ^ e/2 to begin with (recall that /x(Mo) < oo), it 
follows from (2.57) that ||/ — -Pl/||m ^ s. In turn, from this and (2.44) we conclude 
that 

||/--Pl'/||m^£ (2.58) 
for every L' e £(M, 7W, /x) such that L C L' . This finishes the proof of (2.45). 
Next, one notes that the strong convergence in (2.45) implies 

tri,2(M;d/:*)(^) = suptrj,2(M;d^)(Pz,^PL) > suptri2(M;d^)(PLBPi,). (2.59) 

Here we used A ^ 0, the monotonicity of tvL2(^M;dtJi){PLA.PL) with respect to L, 
and (a special case of) the following Lemma 2.8 (with Sn = S = A, Rn =T* = Pl) 
in the first equality in (2.59) and 

J J 
T^iXM, , ^Xm, )m > $^(Xm, , 5Xm, )m (2.60) 

as a consequence of A ;^ S > in the last inequality in (2.59). The non- 
commutative Fatou lemma, [135, Theorem 2.7(d)], and (2.59) then imply B e 
Bi{L'^{M;dfi)), and 

SUptr L2^M;dti){PjBPj) > \\B\\B,(L2(^M;dix)) = ^'^L^(M-4ti){B), (2.61) 

JeN 

where we used S > in the last equality in (2.61), proving (2.36). □ 
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Lemma 2.8. Let p G [l,oo) and assume that R, Rn,T,Tn G B{'H), n G N, satisfy 
s-lim„^oo -Rn = R and s-lim„^oo = T and that S,Sn G Bp{'H), n G N, satisfy 
lim„^oo \\Sn - SWb^ch) = 0- Then lini„^oo ||i?„S'„T* - i?S'T*||Bp(^) = 0. 

This follows, for instance, from [84, Theorem 1], [135, p. 28-29], or [145, Lemma 
6.1.3] with a minor additional effort (taking adjoints, etc.). 

We emphasize that items {i)-{iii) in Lemma 2.7 are not optimal. Indeed, 
one easily verifies that in analogy to (2.4), under the assumptions that A,B£ 
B{L'^{M] dfj.)) and A)^ 0, B )^ 0, one concludes that 

=^ B =^ A if and only if \Bf\ ^ A\f\, f G ^^(M; dn). (2.62) 

On the other hand, the condition < ^|/|, / e L'^{M;dn), in (2.62) remains 

meaningful in the more general situation where B is no longer assumed to satisfy 
B )p Q. In fact, this leads to the following notion of pointwisc domination: 

Definition 2.9. Assume A,B € B {L^ (M ; d/j,)) and A ^ 0. Then A is said to 
pointwise dominate B if 

\Bf\^A\f\, feL\M-dii). (2.63) 

One then has the following improvement of Lemma 2.7 (m), {Hi) due to Dodds 
and Frcmhn [53], Pitt [125], and Simon [135, Theorem 2.13]: 

Theorem 2.10. Assume A,B£ B(L'^{M; dfj,)) , that A)p= 0, and that A pointwise 
dominates B. Then the subsequent assertions hold: 
{i) The following norm hound is valid, 

\\MB{L'2{M-dii)) > \\B\\B{L-^{M;diJ.))- (2.64) 

(m) Suppose in addition that A G B2n{L^{M;d^)) for some n G N. Then also 
B G B2n{L^{M;dij)) and 

\\A\\B2„{L^{M;dp.)) ^ l|S||B2„(i,2(M;dM))- (2.65) 

{Hi) Suppose in addition that A G Bao{L^iM; dfi)) . Then also B G Boo{L'^{M; dfi)) . 

Since obviously, ||B/||m < ||^I/I||m < \\A\\B{LHM;d^.))\\f\\M, f € L^{M;dfi), this 
settles item (i) in Theorem 2.10 (and hence also item (i) in Lemma 2.7). Item 
(a) is proved in [135, Theorem 2.13], and {Hi) is due to [53] and [125] (see also 
[107]). The case of integral operators A and B was first treated in [147, p. 94]. 
We emphasize that Theorem 2.10 (m) is wrong with 2n, n G N, replaced by any 
p < 2 (in particular, it fails in the trace class case p = 1), it is also not true 
for p G (2, (X))\{2(n+ 1)}„£n, due to counterexamples by Peller [122], employing 
Hankel operators (see also the discussion in [135, p. 24, 128-129]). 

Next, we turn to positivity preserving contraction semigroups. First, we recall 

the following well-known result: 

Theorem 2.11 ([128], p. 204, 209). Suppose that H is a semihounded, self-adjoint 

operator in L'^{M; dfi) with Aq = 'mi{a{H)). Then the following assertions {i)-{iii) 

are equivalent: 

{i) e"*^ )pO for allt^O. 

(ii) {H - XIm)^^ > for all A < Aq. 

{Hi) f G dom(|i?|i/2) implies \f \ G dom{\H\^^^) and 

||(F-Ao/m)^/^|/||L< ||(ff-Ao/M)^/V||M- 
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The equivalence of items (i) and (m) follows from the relations 

{H-XIm)~^= dte'^e-^", X < Xq, (2.66) 
Jo 

e-*" = e-*^°s-lim\lM + (t/n){H-XolM)r", t^O. (2.67) 
The next result recalls basic semigroup domination facts: 

Theorem 2.12. Suppose that Hj, j = 1, 2, are semibounded, self-adjoint operators 

in L'^{M;d^) with Xq = iiif {rr { H i), a {H 2))- Moreover, assume that Hj generate 
positivity preserving semigroups, exp(— fifj) )^ 0, j = 1,2. Then the following 
assertions {i)-{iii) are equivalent: 

{i) exp(-iili) )>= exp{-tH2) ^ for all t > 0. 

(ii) [Hi - XIm)-^ > {H2 - XIm)-^ > for all A < Aq. 

(m) For all fj e dom{Hj)+, j = 1,2, (/i,i/2./2)M ^ (i?i/i, ./2)m. 

Suppose in addition that the form domains of Hi and H2 coincide, dom(^\Hi\^^^) = 

dom (|i?2|^^^) • Then items {i)~{iii) are further equivalent to 

(iv) For all fj e dom (|i?i|^/^)_^ = dom (|i?2|^/^) + , j = 1, 2, 

((if2-Ao/M)^/Vi, {H2-XolMY/^f2)j^ > {{Hi-XolMy/yi, {H^-XqImY'^ fl) j^- 

Proof. The equivalence of items (i) and {ii) again follows from the relations (2.66), 
(2.67). 

In order to prove that item (?') imphcs item [Hi] suppose that fj G dom(ifj)_|_, 
j = 1,2, then item {i) implies {fi,er'"\f2)j^j ^ {fi,e^''"'\f2),^ and hence 

{fi,t-^ [hi - e-*^^] f2),, > {fi,t-^ [Im - e-'''^]f2),,. (2.68) 

Letting f 4- in (2.68) yields item {Hi) as e^*^J , j = 1, 2, are strongly differentiable 
with respect to t > 0. 

To show that item {Hi) implies item {ii) one can argue as follows: Let fj G 
L'^{M;dij)+, j = 1,2, then by Theorem 2.11, {H j - XI m)~'^ f j S dom(f/^j)+, A < Aq, 
j = 1,2, and hence by item {Hi), 

{{H^-XIm)-^ fi,H2{H2-XlM)-^ f2) > {Hr{Hr-XlM)-^fi, {H2-XIm)-^ f2) 

(2.69) 

Using H{H - XIm)^^ hi + X{H - Xhi)^^ on either side of (2.69) yields 

(/i, {Hi - A/m)-V2)m ^ {H2 - A/m)- V2)m- (2-70) 
and hence yields item {ii). 

Next suppose that ^ fj G dom (|If2|^^^), j = 1,2, then item (i) once again 
implies 

(/i,i-^[/M-e-*(^-^''^")]/2)^> (/i,t-^[/M-e-*(«-^''^")]/2)^. (2.71) 
Letting t ^ in (2.71) then implies the sesquilinear form version of {Hi), that is, 

{{H2-XolM)'^^fl, {H2-XolM)'^^f2)^ > {{Hr-X<ylM)^l^fl, (i^l - Ao/m)'/'^) m> 

(2.72) 

and hence item {iv). 

Finally, let fj e dom(/fj)+, then fj e dom (|il2|^^^)_,_, j = 1,2, and item {iv) 
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implies 

ifl, {H2 - XolM)f2)M = {{H2 - Ao/m)'/'/i, iH2 - Ao/m)'/'/2)m 

> {{H, - Ao/m)'/Vi, (^^1 - Ao/m)'/'/2)m = ((^^1 - Ao/m)/i,/2)m> 
implying item {Hi). □ 

In the proof of Theorem 2.12 we used the known fact that for any semibounded, 

self-adjoint operator _ff in a complex separable Hilbert space H the domain and 
form domain of H can be characterized in terms of its semigroup e~*^, t ^ 0, by 

dom(iJ) = 1^ e 'H limt"^ II [/-H - e"*^]^!!^ exists finitelyj, (2.74) 
in particular, 

\imt-'[ln-e-'"]f = Hf, / e dom(ff). (2.75) 



44,0 



In addition, 

dom {\H\^/^) = ^ge'H limt'^g, [in - e"*^]^)^ exists finitelyj, (2.76) 
and hence, 

(2.77) 

and by polarization, 
]imt-i(/, [I-H - e-'"]g)^ = {\H\^/' f,sgn{H)\H\^/'g)^, f,ge dom{\H\^/'). 

(2.78) 

The results (2.74)-(2.77) follow from a combination of the spectral theorem for 
self-adjoint operators and Lebesgue's Dominated Convergence Theorem. (In fact, 
(2.74)-(2.78) also extend to unitary groups e~**^, but in this case the proof of the 
analogs of (2.76), (2.77) is more subtle and involves certain results on characteristic 
functions and moments of probability measures, see, e.g., [74, Satz 2.1].) 

Theorem 2.12 is due to [28] and for convenience of the reader we quickly 

reproduced the arguments in [28]. The elementary quadratic form part {iv) in 
Theorem 2.12 is a special case of a general result in connection with m-accretive 
generators of semigroups and their associated forms due to Ouhabaz [117], [119, 
Theorem 2.24], based on the concept of ideals of subspaces of L'^{M; dfi). For 
completeness, we offered the elementary proof of part (iv) based on the observation 
(2.78), which also differs from the arguments used in [137] in the special case of 
symmetric forms. Part (iv) of Theorem 2.12 will subsequently be applied in the 
context of elliptic partial differential operators in divergence form in Section 4. 

Next we turn to additional results on positivity improving semigroups. We recall 
that if no nontrivial. closed subspace of L^(Af ; d/i) is left invariant by T and every 
bounded operator of multiphcation on L'^{M;d^), then this is indicated by the 
statement that {T} U L'^ {M ; dij) acts irreducibly on L'^{M;dij). 

Theorem 2.13 ([128], Thm. XIII.44; [101]). Assume that H is a semibounded, 
self-adjoint operator in L'^{M-, dfi) with Aq = 'mi{(j{H)). Then the following asser- 
tions (i)-(yiii) are equivalent: 

{i) e-*" y for all t > 0. 

{ii) e~*^ >- for some t> 0. 
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{Hi) {e^'^I^^Q U L°°(M;(iy^) acts irreducibly on L"^ {M ; d^) . 

(iv) e~*^ U Z/°°(M; djj) acts irreducibly on L'^{M; djj) for some t> 

[v) {H - Mm)-'^ y for all X < \o. 
{vi) {H — \Im)^^ y for some A < Aq. 

(vii) {{H - XIm)^^} UL°°(M;d^) acts irreducibly on L'^{M;dfi). 

(via) [H — \Im)~^ U L°^{M; dji) acts irreducibly on L?{M; d^) for some A < Aq. 

The following remarkable consequence of Theorem 2.13 will be useful later on. 

Corollary 2.14 ([101]). Suppose that Hj ^ 0, j = 1,2, are semibounded, self- 
adjoint operators in L^{M;dij), and assume that Hj generate positivity preserving 
semigroups, exp{—tHj) )pO, j = 1,2, satisfying 

e-*^i > e-*"^ ^ for all t > 0. (2.79) 

// either e~*^i >- 0, or e~*^^ >- for some t>0, then 

either e'^"^ >- e'*"', or else, e'^"^ = e'*"' for all t > 0. (2.80) 

We also remark that these notions of positivity preserving (resp., improving) 
naturally extend to a two-Hilbert space setting in which one deals with a second 
Hilbert space L'^{Y; dv) with Y d X and Ji = /Lt|y, see, for instance, [28], [101]. This 
is also frequently done in connection with (nondensely defined) quadratic forms (cf., 
e.g., [50, p. 61-62]). 

For subsequent purpose in Section 4 we also recall the following facts on sub- 
Markovian operators: 

Definition 2.15. An operator A e B[L'^{M;dfi)) is called L°° -contractive if 

P/IUoo(M;d^) < ||/|Uoo(M;d^), f € L\M ; d fx)) H L°° {M ■ d tx)) . (2.81) 

Moreover, A is called sub-Markovian if it is positivity preserving, A)^ 0, and L°°- 
contractive. 

Finally, a semigroup {T{t)}t^o C B(L'^{M;diJ,)) is called a contraction semigroup 
on LP{M;d^) for some p > 1, if for all t ^ Q, T{t) extends to a contraction on 
LP{M-dij). 

One verifies that A € B{L-^{M ; d/i)) is i°°-contractivc if and only if it leaves the 
closed, convex set C = {/ G L^{M; dfi) | ]/] ^ 1 /i-a.e.} invariant, that is, AC C C 
(cf. also [115]). Moreover, one readily verifies the following equivalent definition of 
A £ B(^L'^{M;dfi)) being sub-Markovian, namely, 

A is sub-Markovian if and only if / S L'^{M;dn)), < / < 1 

implies Os^Afs^l. ^^'^^"^ 

In fact, (2.82) is used in [66, Sect. 1.4] as the basic definition for a bounded operator 
to be Markovian (an alternative name for sub-Markovian). 

In addition, one notes that if A € B(^L'^{M;dfi)) pointwise dominates the op- 
erator B G B(L^{M;d^)) , and if A is an L°°-contraction, then so is B, that is, 
if 

\Bf\ < A\f\, f e L'^{M;dn) and A is an L°°-contraction, 

(2.83) 

then B is an L -contraction, 
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as the following elementary observation shows: 

ll^/llL°°(M;d^) = \\\Bf\\\L°°{M;dt^) ^ \\A\f\\\L=°(M;dfi) < II |/| ||L°°(M;d/i) 

= ll/l|L~(M;dM), f &L^{M;dfi))nL°-{M;d^i)). 



(2.84) 



Historically, the circle of ideas described in this section originated with Perron 
[124] and Frobenius [63]-[65] in the context of matrices, and was first extended to 
the infinite-dimensional case of compact integral operators by Jentzsch [89] (see also 
[91, p. 183]). The issue, apparently, was revived within the mathematical physics 
community by Glimm and Jaffe [75] around 1970 in the context of nondegenerate 
ground states (i.e., a simple eigenvalue at the bottom of the spectrum) of quantum 
(especially, quantum field theoretic) Hamiltonians. In this connection we refer, for 



instance, to [4], [5], [8], [9], [10], [12], [18], [28], [46], [47, Ch. 7], [51, Ch. 13], [57], 
[58, Sects. 8, 10], [59], [62], [67], [76, Sect. 3.3], [77], [81], [87], [88], [101], [102], 
[105], [106], [114], [115], [117], [119, Chs. 2, 3], [123], [126], [128, Set. XIII.12], [129], 
[130], [131], [132], [134], [143, Sect. 10.5], [146, Sect. 3.3] and the references cited 
therein for the basics and some of the applications of this subject. 



3. Robin- Type Boundary Conditions for Matrix- Valued Divergence 



In this section we develop the basics for divergence form elliptic partial differen- 
tial operators with (nonlocal) Robin-type boundary conditions in bounded Lipschitz 
domains. In fact, we will go a step further in this section and develop the theory 
in the vector-valued case as this is certainly of interest in its own right. Thus, we 
will focus on m X m, m e N, matrix-valued differential expressions L which act as 



For basic facts on Sobolev spaces and Dirichlet and Neumann trace operators, as 
well as the choice of notation used below, we refer to Appendix A. For the basics 
on sesquilinear forms and operators associated with them we refer to Appendix 
B. Moreover, for the definition of Lipschitz domains ft and the associated Sobolev 
spaces on fl, i?*(ri), and its boundary dfl, H^{dfl), we refer to Appendix A. 

In the remainder of this section we make the following assumption: 

Hypothesis 3.1. Let n e N, n > 2, and assume that fl cW^ is a bounded Lipschitz 
domain. 

For simplicity of notation we will denote the identity operators in L'^{Q;d"x) 
and L^(90; rf"~^w) by In and Ign, respectively. Also, in the sequel, the sesquilinear 
form 



( • , • = H3(aQ)-( • , • : H'{dfir x H-^dQr ^ C, s e [0, 1], 



(antilinear in the first, linear in the second factor), will denote the duality pairing 
between H^{dCl) and 



Form Elliptic Partial Differential Operators 




(3.1) 



(3.2) 



H-'idn)"" = {H'{dnD*, s e [o, i]. 



(3.3) 
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such that 

JcKi (3.4) 

/ e H'idn)"', g e L^idn- d^-^w)" ^ H-'idn)"", s e [o, i], 

where d"'~^oj stands for the surface measure on d^l. 

Below we introduce the class of multipliers on generic Banach spaces (though 
we are primarily concerned with Sobolev spaces; for an authoritarian treatment of 
this topic the interested reader is referred to [112]). 

Definition 3.2. Let (X, || • \\x), {Y, \\ • \\y) be two given Banach spaces of distribu- 
tions in a domain f2 C M" with the property that C°°{Q) is dense in both of these 
spaces. In this context, define A4{X Y), the space multipliers from X to Y , as 
follows: If Ma denotes the operator of pointwise multiplication on C°°{Q) by the 
function a, then a € J^{X ^ Y) indicates that Ma may be extended to an operator 
in B{X,Y). Whenever a G M{X Y), we set 

\\a\\M{X^Y) ■■= \\Ma\\B{X,Y)- (3.5) 

Finally, we abbreviate M{X) := M{X ^ X). 

Next, we wish to describe a weak version of the normal trace operator associated 
with an m X m, to G N, second-order system in divergence form (3.1), considered 
in a domain O C M". To set the stage, assume Hypothesis 2.1 and suppose that 
some s e (0, 1) has been fixed with the property that the tensor coefficient 

^ {(^j,k) Kj,fc<« (3.6) 

l^a,/3^m 

satisfies 

Ae MiH-'-^/^n)""). (3.7) 
One observes that the inclusion 

t: H'^'in) ^ (H^'in))*, So > -1/2, r > 1/2, (3.8) 
is well-defined and bounded. We then introduce the weak Neumann trace operator 

^N- {uG H'+^/'^{n)"'\Lue H'^in)""} ^ H'-\dnr, so > -1/2, (3.9) 

as follows: Given u e i?"+i/2(fi)"' with Lu e i?"«(fi)'" for some sq > -1/2 and 
s e (0, 1), we set (with t as in (3.8) for r := 3/2 - s > 1/2) 

{(P,1NU)^^1 = jjl/2-s(Q)m(£'$, A£)u)(jjl/2-.(f2)™). 

- jj3/2-»(Q)m($,f.(iu))(jj3/2-3(Q)m)., (3.10) 

for all (t> e ffi-«(90)'" and $ e if3/2-s(n)™ such that jd^ = (j). Above, we used 

Du := [djUc)l4^a<^m (3-11) 

to denote the Jacobian matrix of the vector- valued function u = (wa;)i^a^Tn) and 
we employed the convention that, in general, 

n rn 

■= (E E <f ) for all C = {C',U,^m € C"x™. (3.12) 

fc=l 3=1 l^j^n l^fc^n 
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In addition, we will later employ the eonvention 

n m 

{V,AC):= Y: E «;.fCf<, (3.13) 

j,k=l a, 13=1 

for all r? = e C = (Cf )K/3^m G C"><™. 

Returning to the mainstream discussion, we note that the first pairing on the 
right-hand side of (3.10) is meaningful granted (3.7), since 

(/f^/2-"(0)")* = H'-'^/^{Q)'". (3.14) 

Moreover, the definition (3.10) is independent of the particular extension $ of (f>. 
Indeed, if $' e _ff"^/^~®(i7)™ is another vector-valued function such that jd^' = 4> 
then 4- :=$-$' e H^^'^-%Q)"' satisfies 713* = 0. Hence * G HQ^^^'{n)"^ and, 
as such, there exists a sequence {tl>j}jeN C CQ°(f2)'" with the property that 

V'j -s- * in H^/'^-^{Q)"', as j -s- 00. (3.15) 

Consequently, given that Dipj in H^/'^~^{Vl)'^ as j 00, we may compute 

/J-a/^-s (Q)"> (-(iw)) (ij3/2-s(Q)m). 

= ^.1™^ /J-3/2-'*(n)"' (V'j) ^'{LU)) (^H3/2-s 

J—KX> 

= ^Ifin x,(n)m(£)V'j,-4Du)(p(n)m)/ 

= jji/2-a(n)-(^*,>l£'M)(//i/2-=(n)-)*, (3-16) 
with X>(il) denoting the space of test functions (i.e., CQ°(il)) equipped with the 
usual inductive limit topology (so that, in particular, V^il) =(C^(r2)) is the space 
of distributions in O). Now (3.16) gives that jji/2-s(^Q-jm.{D'^ , ADu) — 
//3/2-s(Q)™ (^', t(LM))(^3/2-3(Q)m). = which ultimately shows that 

Hi/2-,(n)™(D<I>, AI?M)(^i/2-«(n)™). - //3/2-»(n)™(^:'-(-^w))(//3/2-«(n)™)* (3-17) 

Prom this, the desired conclusion (pertaining to the unambiguity of defining 7jvU 
as in (3.10)) follows. 

One recalls that has a linear right-inverse, that is, there exists a linear and 
bounded operator 

£ : H'{dfl) ^ H'+'^/'^{fl), 0<s<l, (3.18) 

which is universal (in the sense that it does not depend on s G (0, 1)) and satisfies 

7D(^0) = (?i, (j>€H'{dn), 0<s<l. (3.19) 

Given an arbitrary (j) G H^{dO,) and choosing $ := f e H^+'^/'^{Q) in (2.9) then 
allows us to estimate 

|l7Aru||H»-i(ao)- s; C'(P||A^(ffa-i/2(n)™)||?i||ff,,+i/2(n)™ + ||Lu||^3o(n)™) (3-20) 

for every u in the domain of jn- This proves that the operator in (3.10) is 
well-defined, linear, and bounded. 
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It is instructive to point out that, in the case when s = 1/2, condition (3.7) 
reduces precisely to 

AeL°°(f2;rf"a;)'"^™. (3.21) 

In particular, 

^N:{u€H\nnLu€H'"{nr}^H-^/^{dnr, So > -1/2, (3.22) 
is well-defined, linear, and bounded, whenever (3.21) holds. 

Hypothesis 3.3. Assume Hypothesis 3.1, suppose that 6 > is a given number, 
consider m G N, and assume thatee B{H^''^{dn)'^ ,H-^/'^{d^)'^) is a self-adjoint 
operator {in the sense discussed in (B.7)) which can he written as 

e = ew + e(2' + e(3\ (3.23) 

where Q^^\ j = 1, 2, 3, have the following properties: There exists a closed sesquilin- 
earform qf^ in 1.2(90; d^-^w)™, with domain H^'^idCl)"' x H^/'^{dCt)"', which is 
bounded from below by cqq G R {hence, q^j'^j is symmetric) such that if Qg'^ ^ 
Cdnldn denotes the self-adjoint operator in L'^{dQ.\d"'^'^uj)"^ uniquely associated 
with (c/. (B.27)), then 9^^^ = ©qq, the extension of to an operator in 
B(ijV2(an)™,if-i/2(5f2)m^ discussed in (B.26) and (B.32)). In addition, 

6(2) e B^{H^'^{d^)'^ ,H-^'^{dQ.)"'), (3.24) 
whereas G^^) g B{H^/'^{dn),H-^/'^{dn)) satisfies 

ll®''^^llB(ffi/2(an)",(_f/-i/2)™(ao)) < (3.25) 

We record the following useful result involving the Dirichlet trace operator 713. 

Lemma 3.4. Assume Hypothesis 3.1 and fix m & N. Then for every £ > there 
exists a /3(e) > {with /3(e) = 0(l/e)) such that 

(3.26) 

Proof. The case when m = 1 (corresponding to the situation when u is scalar- 
valued and hence Du = Vu) has been treated in [69] . The vector- valued case then 
follows from this by summing up such scalar estimates. □ 

Lemma 3.4 is a key ingredient in proving the i7^(ri)-coercivity of the sesquilinear 
form 0e,n,L( • , • ) Theorem 3.6 below. Before stating it, we recall our conventions 
(3.12) and (3.13) and introduce the following hypothesis. 

Hypothesis 3.5. Assume Hypothesis 3.1 and suppose that L is an m x m, m (1 M, 
second-order system in divergence form expressed as in (3.1) for a tensor coefficient 
as in (3.6). In addition, assume that (3.21) holds, that is, A € L°°(f2; d"x)'"'^'", 
and that the tensor coefficient A satisfies the strong Legendre ellipticity condition 
{for some ao > 0) 

(n m \ 

E E «",fCfg UaoKr, (3.27) 
j,k=l a,/3=l ^ 
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for every ( — ((j) i^r^n G C"^™. Finally, suppose A satisfies the symmetry 

condition 

<f = <'J' j,ke{l,...,n}, a,pe{l,...,m}. (3.28) 

The following result extends work from [69] carried out for the special case of 
the Laplacian —A. 

Theorem 3.6. Assume Hypothesis 3.3, where the number S > is taken to be 
sufficiently small relative to the Lipschitz character of fl, more precisely, suppose 
that < 5 < i/i/2(9n))' addition, assume Hypothesis 3.5 and 

consider the sesquilinear form Qe,n{ • , • ) defined on iJ^(O)'" x 77^(0)™ by 

Qe,n{u,v):= f d"x{{Du){x),A{x){Dv){x)) + {jDU,ejDv)[":l ^ ^ 
Jn I (3.29) 

e iJi(n)'". 

Then there exists k > 1 with the property that the form 

i3e,n,K(M,?^) := iJe,n(w,w) + K(w,w)L2(n;<i"a;)™, u,v & H'^{n)"\ (3.30) 

is H^ (Q)"^ -coercive. As a consequence, the form O0,n( • , • ) in (3.29) is symmetric, 
H^{Q)"^ -bounded, bounded from below, and closed in L'^{fl;d"'x)"^ . 

In the proof of Theorem 3.6, the following abstract functional analytic result is 
going to be useful. 

Lemma 3.7 ([69]). Let V be a reflexive Banach space, W a Banach space, assume 
that K e Boo(V, V*), and that T e B{V,W) is one-to-one. Then for every e > Q 
there exists Ce > such that 

\v{u,Ku)v'\<e\\u\\l + C4Tu\\l^, u&V. (3.31) 

We are now ready to present the 

Proof of Theorem 3.6. We shall show that k; > can be chosen large enough so 

that 

\Ml.(a)m ^\ j^d^x\{Du){x)\^ + ^ jj^x\u{x)\'' + {mQ^^hDuf;^l 

u€ H^in)"", j = 1,2,3, (3.32) 

where G^-'^ j ~ 1,2,3, are as introduced in Hypothesis 3.3. Summing up these 
three inequalities then proves that the form (3.30) is indeed iJ^(fl)'"-coercive. To 
this end, we assume first j = 1 and recall that there exists cso S M such that 

(7D^i,e(i)7i3«)!7i > ce„||7DM||i2(aj,^,„-.,)„, u G H\nr. (3.33) 
Thus, in this case, it suffices to show that 

1 

< ^ / d"a;|(DM)(x)|2 + J /" rf"a;|u(a;)|2, ueH^n)"", (3.34) 
or, equivalently, that 



max {-ceo > 0} ll7r>wlli2(ao.d„-i^)„ + ^\\u\\hi 



max{-C0o , 0} ||7r>w|||2(aQ.rf„-i^^)m 



22 F. GESZTESY, M. MITREA, AND R. NICHOLS 

< ^ / d"x|(£)u)(a;)|2 + ^^^ / d»;c|u(:c)|2, uGH^n^, (3.35) 
6 Q Jn 

with the usual convention, 

MUn^r. = l|i^«ll^,.,„,„^,. + Ml.,^.„, u e H\nr. (3.36) 

Now, the fact that there exists k > for which (3.35) holds follows directly from 
Lemma 3.4. 

Next, one observes that in the case where j = 2, 3, the estimate (3.32) is implied 

by 

\{^DuM'hDu)^;^l\ d-x\{Du){x)\'' + '^-^ j d^x\u{x)\\ u€H\nr, 

(3.37) 

or, equivalently, by 

|<7i.n,e«7Dw)S7i| ^ J||«fHi(n)^ + '^MUn;d-^)-' « ^ ^^(O)™. (3.38) 

When j = 2, in which case O'^) e B^{H\n)"^, {H\n)"^)*), we invoke Lemma 
3.7 with 

V := H^n)"", W := L\n, d"a;)'", (3.39) 

and, with 7D G B{H^{n)"',H^/'^{dn)"') denoting the Dirichlet trace (acting com- 
ponentwise) , 

K ■.= jh&2lD e B^{H\Qr, {H\nry), T := ^ : H^n^ L\Q,rxr, 

(3.40) 

the inclusion operator. Then, with e = 1/6 and k := 3Ci/6 + Ij the estimate (3.31) 
yields (3.38) for j = 2. 

Finally, consider (3.38) in the case where j = 3 and note that by hypothesis, 

\{lDU,e^^'^-/Du)':^J^\ ^ ||0^'^||e(Hi/2(ao)-.H-i/2(ao)™)ll7Du||^i/2(an)". 

^ShD\\l^HHn)^M^/Hon)'^)\MHHn)-, u & H\nr . (3.41) 
Thus (3.38) also holds for j = 3 if 

< ^ ^ ^ll7DllB(Hi(n),Hi/^(8Q)) and k>1. (3.42) 

This completes the justification of the estimate (3.32). In turn, this further implies, 
with the help of the strong Lcgcndre ellipticity condition, that 

^II^Uhho)- ^ ^ l^d-x({Du}{xj,A{x){Du){x)) 

+ I ^ \u{x)\^ + {jDU, e(^hDu)^;^l (3.43) 
u £ H\n)"\ j = 1,2,3, 

where A > is as in (3.27). The estimate (3.43) establishes the claim about 
the sesquilinear form in (3.30). Moreover, the symmetry of the sesquilinear form 
£}e,n,L(' , • ) from (3.29) is a direct consequence of (3.28). Finally, the remaining 
claims in the statement of Theorem 3.6 are implicit in what we have proved so far, 
and this finishes the proof of Theorem 3.6. □ 
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Next, wc turn to a discussion of the realization of an m x m system L as a 
self-adjoint operator on L^(il; when equipped with certain nonlocal Robin 

boundary conditions in a bounded Lipschitz subdomain of M". Below, 7jv denotes 
the weak Neumann trace operator discussed in (3.9)-(3.22). 

Theorem 3.8. Assume Hypothesis 3.3, where the number 6 > is taken to be 
sufficiently small as in Theorem 3.6. In addition, assume Hypothesis 3.5. Then 
Le.n, the L'^ -realization of L equipped with a nonlocal Robin boundary condition in 

L^{n;d'^x)"\ defined by 

Le,n = L, 

dom(ie,n) = {w e 77^(17)™ | Lu G L^(f7; d"a;)'", (3.44) 

(7jv + Qid)u = m ^--1/2(90)™}, 

is self-adjoint and bounded from below. Moreover, 

dom{\Le,n\^/^)=H\nr, (3.45) 

and Lq^q,, has purely discrete spectrum bounded from below, in particular, 

(Tess{Le,n) = 0. (3.46) 

Finally, Lq q is the operator uniquely associated with the sesquilinear form Oe,n 
in Theorem 3.6. 

Proof Denote by £2e,a( ■ ) ■ ) the sesquilinear form introduced in (3.29). From 
Theorem 3.6, we know that Oe,n is symmetric, iJ^(fi)™-boundcd, bounded from 
below, as well as densely defined and closed in L"^ {fl; d" x)"^ x (fl; d" x)"^ . Thus, 
if as in (B.34), we now introduce the operator Le,a in L'^{^; d^x)"^ by 



there exists some Wy G L^(0; d^x)™ such that 



dom(i©,n) = jt; e H^{fiy 

J drx(p^,A{Dv)) + {jDW,e'jDv)^^^ = J d'^xw-Wy for all w e ifi(17)™|, 

L0,nu = Wu, u€ dom(Z/e,a), (3.47) 

it follows from (B.20)-(B.35) (cf., in particular (B.27)) that L&^q is self-adjoint and 
bounded from below in 1/^(12; d^x)"^ and that (3.45) holds. Next we recall that 

iJi(f])™ = {m e i?^(n)™ I 7i3U = on dn), (3.48) 

where, as usual, the Dirichlet trace operator 70 acts componentwise. Taking v S 
C^(f7)™ i?o(f^)" ^ i?Hri)™, one concludes that if u e dom(Le,n) then 

/ d"xv-Wu= I d"xv-Lu for all ^; G (17), hence Wu = Lu in (^'(f])™)'. 
Jn Jq 

(3.49) 

Going further, suppose that u G dom(Le,a) and v G H^{fl)"^. We recall that 
7r,: iJi(n)™ ijV2(gn)m boundedly and compute 

/ d"x(D^,A{Du)) = [ d'^xv- Lu+{-yDV,^Nu)'{?^ 
Jq Jq 

= j d'^xv ■ Wu + {'Jdv, (7jv + 67r()u)[™2 - (idv, Q'Jdu) 



(m) 
1/2 



24 F. GESZTESY, M. MITREA, AND R. NICHOLS 

= ^ (Tx (D^, A{Du)) + {jDV, (7iv + Qjd)u)["^^, (3.50) 

where we used the second hne in (3.47). Hence, 

(iDV, (7JV + Q1d)u)'-;J^ = 0. (3.51) 

Since v e iJ^^^)™ is arbitrary, and the map 7d : -ff^^)" ^ i?^/^(aO)'" is actually 
onto, one concludes that 

(7JV + Q7d)u = in H-^/^idn)"". (3.52) 

Thus, 

dom(L0,n) C {i; e H^{n)"' \ Lv e {n; cP x)"^ , 

{jN + ejD)v = in H-'^/^idny''}. 

Next, assume that u e {v e H\n)"^ \ Lv S L^^f^; d"a;)™, (7^ + &-/d)v = O}, 
«; e H^iCl)"", and let w„ = Lu e ^^(n; rf"a;)™. Then, 

(Pxw ■ w,, = / (Pxw ■ Lu 



(3.53) 



/ 



= ^ rx (Dw, A{Du)) - {^DW, ^Nu)^^l 



(3.55) 



crx{Dw,A{Du)) + (7z,w;,e7i5u);"2^ (3.54) 

Thus, applying (3.47), one concludes that u E dom(— Ae,j2) and hence 

dom(Le,n) ^ {v G ^^(O)™ | Av e L'^ {Vl- x)"^ 

{in + Qid)v = in H-^/^{d^ir]. 

Finally, the last claim in the statement of Theorem 3.8 follows from the fact that 
if i(f2)™ embeds compactly into L'^{Vt- d"a;)™ (cf., e.g., [55, Theorem V.4.17]). □ 

Remark 3.9. Wc emphasize the explicit form of the domain of the; operator Le,o 
displayed in (3.44) in terms of boundary trace operators 71) and ^n- This particular 
feature of an explicit domain rather than an operator defined via the underlying 
sesquilinear form and the First Representation Theorem is one of the reasons for our 
choice of Lipschitz domains fJ; it also dictates our conditions on 6 in Hypothesis 
3.3. 

In the special case 6 = 0, that is, in the case of Neumann boundary conditions, 
we will also use the notation 

Qjv,n(-, •) = ^o,q(-, •)> LN,n-=Lo,n- (3.56) 
When specialized to the case m = 1 and L = —A, Theorem 3.8 yields a family 
of self-adjoint Laplace operators — Ae^n in L-^{i};d"x) indexed by the boundary 
operator 9, which we shall refer to as nonlocal Robin Laplacians. More specifically, 
we obtain the following result first proved in [69] . 

Corollary 3.10. Assume Hypothesis 3.3 (with m = 1), inhere the number S > is 
taken to be sufficiently small as in Theorem 3.6. Then —Aq^q, the nonlocal Robin 
Laplacian in L^(0;c?"a;) defined by 

- A@n = -A, 
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dom(-Ae,n) = {u e H\n} \ Au e L^{n; d"x), (3.57) 
{iN + Qjd)u = in 
is self-adjoint and bounded from below. Moreover, 

dom(|- Ae,nr/') =i?^(fl), (3.58) 
and — Ae^n, has purely discrete spectrum bounded from below, in particular, 

(Tess(-Ae,n) = 0. (3.59) 
Finally, — Ae^Q is the operator uniquely associated with the sesquilinear form 

qe,o(M, v) d"x (Vu)(a;) • (Vw)(x) + {-fou, e-fDv)^^, u,v e H\n). (3.60) 

In the special case 6 = 0, that is, in the case of the Neumann Laplacian, we will 
also use the notation 

qAr,f2(-, •) = qo,J2(-, •); -Aat^h := -Ao,o. (3.61) 

Next, we briefly comment on the usual case of a local Robin boundary condition, 
that is, the case where 9 is the operator of multiplication M$ by a function 9 defined 
on d^l: 

Lemma 3.11 ([69]). Assume Hypothesis 2.1 and suppose that 8 = Mo, the op- 
erator of multiplication with a measurable function : dfl K. Suppose that 
9 e LP{dQ] d^'-'^Lo), where 

p = n—l if n > 2, andpG{l,oo] if n = 2. (3.62) 

Then 

&eBoc.{H^^^{dn),H-^/'^{dQ)) (3.63) 
is a self-adjoint operator which satisfies 

ll©llB(i/i/2(ao),_H-i/2(9n)) ^ C'll^'l|Lp(af2;d"-iw)7 (3.64) 

for some finite constant C = C{^,n,p) > 0. In particular, the present situation 
9 = Mo subordinates to the case 9^^^ described in (3.24). 

It is worth noting that the conditions isolated in Hypothesis 3.3 permit one to go 
beyond the assumption 9 G L°°(dil; d'^^^oj) one frequently finds in the literature in 
connection with local Robin boundary conditions and hence naturally lead to the 
iP-conditions in Lemma 3.11. 

In the case 8 = Mg described in Lemma 3.11, the underlying sesquilinear form 
and operator will bo denoted by and 

The L^-realization of L equipped with a Dirichlet boundary condition, Lu^q, in 
L^(f7; d^x)"* formally corresponds to 8 = oo and so we isolate it in the next result. 

Theorem 3.12. Assume Hypothesis 3.5. Then L^ ji, the version of L equipped 

with a Dirichlet boundary condition in {Q.; d"' x)"^ , defined by 

Lo.n = L, 

dom{LD,n) = {u e H^i^)"" \ Lu G L^{n- d"x)"\ jdu = m H^/^idn)""} 

(3.65) 

= {u€ H^iftr I Lu e L^iQ; d^x)""}, 
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is self-adjoint and strictly positive. Moreover, 

dom((Lz?.n)'/') =i^o'(^^^. (3-66) 
and since SI is open and bounded, Ljj q has purely discrete spectrum contained in 
(0,00), in particular, 

aessiLD,a) = 0. (3.67) 
Finally, L^.n is the operator uniquely associated with the sesquilinear form 

QD.n{u,v) := [ d"x{{Du){x),A{x){Dv){x)), w, ?; e Fo^(f7)™. (3.68) 
Jn 

This is largely proved as in thc^ case of Theorem 3.8, and hence we omit the 
argument. Here we only wish to note that (3.67) follows from (3.66) since iJQ(f2)'" 
embeds compactly into L^(0; fi"a;)™; the latter fact actually holds for arbitrary 
open, bounded sets fl C K." (see, e.g., [55, Theorem V.4.18]). 

By specializing Theorem 3.12 to the situation when m = 1 and L = — A yields 
the following corollary (cf. also [68], [71] and [69] for related results): 

Corollary 3.13. Assume Hypothesis 3.1. Then —An,n, the Dirichlet Laplacian 

in L'^i fl: d'."-x) defined by 

- Adsi = -A, 

dom(-Ai3,n) = {m e H\n) \ Au e L^{n; d'^x), -^du = m H'^'^{dn)) (3.69) 

= {u e Hl{n) I Am e L^{n-d'^x)], 

is self-adjoint and strictly positive. Moreover, 

dom ((-Ai,,a)i/2) = HliSl), (3.70) 

and since Q. is open and bounded, — A^^n has purely discrete spectrum contained in 
(0,00), in particular, 

aess(-AD,n) = 0. (3.71) 
Finally, —Ad,q is the operator uniquely associated with the sesquilinear form 

c\dMu,v) ■■= I d"x(Vw)(x) • {Vv){x), u,v€Hl{n). (3.72) 
Jn 

4. Gaussian Heat Kernel Bounds for Divergence Form Elliptic PDOs 
WITH Robin-Type Boundary Conditions 

In this section we apply the abstract results of Section 2 to the concrete cases 
involving elliptic partial differential operators in divergence form in L^(0;d"a;) on 
bounded, connected Lipschitz domains Q. with Robin-type boundary conditions on 
dSl studied in Section 3. 

Throughout this section we consider the scalar case m = 1 and assume Hypoth- 
esis 3.3 whenever a (nonlocal) Robin boundary condition is involved. Moreover, we 
introduce the following assumption: 

Hypothesis 4.1. Letn eN, 2. 

(i) Assume that O c K" is a bounded Lipschitz domain {cf. Appendix A). 

(ii) Suppose that the matrix A{-) is Lebesgue measurable and real symmetric a.e. 
on Q.. In addition, given Q < < a\ < 00, assume that A satisfies the uniform 
ellipticity conditions 

aoln ^ -^{x) ^ ailn for a.e. x G Q. (4.1) 
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Here In represents the identity matrix in C". 

Given the basic setup for divergence form elhptic partial differential operators 
-^e,n, with nonocal Robin boundary conditions developed in Section 3, we can 
now turn to (Gaussian) heat kernel and Green's function bounds for Lq^q (and 
hence for the special case — Ae^n), and subsequently also for the corresponding 
Schrodinger-type operators. 

We will use the following heat kernel notation (for t > 0, a.e. x,y E fl) 

and similarly for Green's functions (for z € C\K, a.e. a:,y G fi), 

Ge,n(2,a;,y) = (ie,n - z^n)~\a;, y), GN,n{z,x,y) = {LN,n - zla)~^{x,y), 

GD,a{^,x,y) = {Ld,q- zln)~^{x,y), x^y. (4.3) 

Next, we need one more preparatory result which identifies i_D,n as the strong 
resolvent limit of Le.fj as ^ oo. For this purpose we recall an ordering < on the set 
of nonnegative quadratic forms in a complex, separable Hilbert space "H as follows. 
If ti and t2 are two non- negative quadratic forms in 'H, then 

t2 <ti if and only if dom(ti) C dom(t2) and iii^u) ^ ti(M), u G dom(ii). (4.4) 

Monotonically increasing sequences of non-negative quadratic forms have limits. 
More precisely, one has the following fact (cf. [29, Lemma 5.2.13], [94, Ch. VIII, 
Theorem 3.13a], [132], [133]): 

Lemma 4.2. Suppose that {tn}neN denotes a monotonically increasing sequence 
{with respect to the ordering < for quadratic forms introduced in (4.4)) of non- 
negative, closed quadratic form„s in a complex, separable Hilbert space H. Then the 
following items (i) and (ii) hold, 
(i) i defined by 

i{u) = supin{u), 



u e dom(l) = iuGl-L u e Pi dom(t„), supt„(u) < oo >, 
I LI nen J 



(4.5) 



is a non-negative, closed quadratic form in T-L, and one has 

lim t„(u, u) = t(?i, w), M,wGdom(t), (4.6) 

n— >oo 

where t,i(-, •), n G N {resp., t(-, •)) is the sesquilinear form defined from the quadratic 
form tn, n gN {resp., t) via polarization. 

{ii) Suppose in addition that i {and therefore each t„, n G N) is densely defined. 
If Tn, n G N, {resp., T) denotes the unique densely defined, non-negative, self- 
adjoint operator associated to the sesquilinear form t„(-, •) {resp., t(-, •)) by the First 
Representation Theorem, then 

s-lim (T„ - XI-h)-^ = {T- XIn)-\ A < 0. (4.7) 

n— >oo 

In particular, 

s-lim e"*"^" = e"*"^, t > 0. (4.8) 
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Wc note that item (i) and (4.7) are taken from [29, Lemma 5.2.13] and [94, Ch. VIII, 
Theorem 3.13a]. By analyticity, (4.7) extends to all A e C\IR (cf., e.g., [143, 
Theorem 9.15]), 

s-\im{Tn-zlH)~^ = {T-zIh)-\ ^eC\[0,oo). (4.9) 

n— ^oo 

Since T„ , n <E N, and T are non-negative, choosing z = i{= \/— 1) in (4.9) and 
applying [143, Theorem 9.18] yields (4.8). 

The following result is mentioned in [101] in the context of Dirichlet Laplacians 
and we thank Derek Robinson for helpful discussions concerning its proof: 

Lemma 4.3. Assume Hypothesis 3.3 (with m = 1) and denote by L^,q the operator 
Le,n ii^ the special case where O denotes the operator of multiplication, M^, in 
L'^{dfl;d"~^ui) with the positive constant ?9 > 0. Then 

s-lime"*'^''''^ = e"*^"-'', t^O. (4.10) 

Proof. Let {'&n}?^=i denote an arbitrary sequence of positive real numbers satisfying 
< i?„ < iJn+i, n e N, and lim ■!?„ = oo. (4-11) 

n—>oo 

The inequalities in (4.11) imply the sequence {Qs„,q{-, •)}^=i of non-negative, 
closed quadratic forms (cf. (3.29) with O replaced by M^^) is monotonically in- 
creasing with respect to the quadratic form ordering < introduced above in (4.4). 
In this simple case, the quadratic forms have a common domain (viz., H^{fl)); thus, 
monotonicity simply amounts to 

ntf„,n(w, u) ^ 0,?„+i,n(u, u), u e H\n), n e N. (4.12) 

By Lemma 4.2, 

Q(x),n('«, w) = sup0tf„,n(u,u), u e dom(£2oo,n(-, •))> 

neN 

dom(Ooo,a(-,-)) (4-13) 
= \ u € {CI; d" x) u€ n dom(i3^^,n (•,•)), supi!3^„,n('U, u) < cxd >, 

defines a non-negative, closed quadratic form in L^( fl; d^x). Wc claim the limiting 
quadratic form 0oo,o(', •) defined by (4.13) coincides with Q.d,o,{', the quadratic 
form of the Dirichlet Laplacian LD,n in L'^{Cl;d^x). In order to prove this, one 
needs to verify the following two conditions, 

Aom{Q^^n{-.-)) = Hl{n), (4.14) 
O.oo,n{u, u) = O.D,n{u, u), u G Hl{H). (4.15) 

To verify (4.14) one observes that u G dom(Ooo,n(-, •)) if ^^'^ o^^y u G H^{Q) 
and 



sup 

neN 



/ d^'x {{Vu){x),A{x){Vu){x)) + K{JDU,^DU) 

Jn ' 



< 00. (4.16) 



1/2 

Clearly, (4.16) is fulfilled if and only if 
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that is, if and only if 7_dm = in H^/^{Q). Therefore, u G doni(Ooo,n(-, •)) 
and only if m e H^{Cl) and 7dU = in and (4.14) follows. Finally, one 

computes 

noo,n(w,w) = supOtf„,n('W, w) = / d^-x {{Vu){x),A{x){Vu){x)) = Q.d,u{u,u), 
neN Jn 

uGH^{n), (4.18) 

implying (4.15). By Lemma 4.2, and (4.8) in particular, with r„ = L^„,q, n gN, 
and T = Ld,q, 

s-lime"*^''"'" = e"*^°'", t ^ 0. (4.19) 

Since {??n}J^Li satisfying (4.11) was arbitrary, (4.10) follows. □ 

Now we arc ready to formulate the first principal result of this section. 

Theorem 4.4. Assum,e Hypothesis 4.1. suppose that Oj, j = 1,2, satisfy the 
assumptions introduced in Hypothesis 3.3, and denote by ie^.Ji the operaior's in 
(3.44) uniquely associated with the sesquilinear forms l}ej,o( ' j ' ); J = Ij^; defined 
on H^{ft) X H^{fl) according to (3.29). Suppose, in addition, that 

<7D|w|,e,7z3|u|)^/2 ^ {jDU,&jjDu)^^^, u&H\n), J = 1,2. (4.20) 
Then, assuming ^ Oi ^ 82, one has the positivity preserving relations 

=^ e-*^°'" =^ e-*^«2." =<; Q-tLo^.n ^ e"*^"'", t > 0, (4.21) 
or equivalently, 

=^ (Ld,q + A/n)-^ =^ (Le2,a + A/q)"^ =^ {Le„a + A/q)"^ 

1 (4-22) 
=^ (iiv.n + A7o)-\ A>0. 

In addition, all semigroups appearing in (4.21) lie in the trace class, 

e-tLn.n^ g-*ie,,f!^ g-tiiv,n g Bi{L^ i^; d"" x)) , j = 1,2, t > 0. (4.23) 

In particular, one has the Gaussian heat kernel bounds (for t > 0, a.e. x,y € Cl), 

^ KD,n{t,x,y) ^ K0^,i2{t,x,y) < K^^^i^it, x,y) ^ K]y,ii{t,x,y) 

/o \ r n . (4-24) 

«;C„,„„,omax(t-"/2,l)exp{-|a;-y|V[4(l+7)aii]}, 7^(0,1). 

anrf Green's function bounds (for A > 0, a.e. x,y d fl), 

< GD,a(A,x,j/) ^ Ge2,nit,x,y) ^ G'ei,n(A, a:, y) ^ GN,a{\x,y) 

^ fCao,ai,A,n,nk-yp-", n > 3, (4.25) 

^ \C'ao,ai,A,a|ln(l + |a;-2/|-^)|, n = 2, 

Proof. By [50, Theorems 1.3.5, 13.9], Xd^q and ijv,n satisfy the Beurling-Deny con- 
dition {Hi) of Theorem 2.11, and hence by Theorem 2.11 (i), e"*^" " and e~*^"'", 
t ^ 0, are positivity preserving. In addition, both have nonnegative integral kernels 
(cf. also Theorem 2.3) and are known to satisfy the bounds (C.3) and (C.4). 

As the relations (4.21) and (4.22) are equivalent by Theorem 2.12, it suffices 
to focus on (4.21). But then, =^ e'^^'^^.^^ ^ g-tiei.si ^ g-tLiv.r;^ t ^ 0, is 
immediate upon combining Theorem 2.12 (iv), (3.29), and (4.20). Moreover, using 
Lemma 4.3, (4.21) follows. Thus, an application of (2.15) implies the inequalities 
(4.24) and (4.25), except, the very last in either one of them. The inequality 
(C.4) for KN^Q{t,x,y) then completes the proof of (4.24). Similarly, the inequality 



30 



F. GESZTESY, M. MITREA, AND R. NICHOLS 



Gjv,n(A,a;,y) ^ Can, ai, a, o,n|a; - X > 0, x,y e il, x ^ y, can be found, for 

instance, in [54, Lemma 3.2] (see also [35]) for n ^ 3. As wo were not able to find 
the case n = 2 in the literature, we provide a short proof in Appendix C. 

The inequalities (4.24) then follow from Theorem 2.3, (4.21), and (4.22), as soon 
as we establish that semigroups of the type with 8 satisfying the conditions 

of Qj are, in fact, integral operators. This follows from combining Lemma 2.7 {iii) 
and (4.21) as the Neumann heat kernel bound (C.4) yields the trace class property 
g-tiiv,f2 g Bi(L^(f2;d"a;)), t> Q (this trace class property of course also applies to 
e-*^° ", t > 0). 

Finally, using (2.66) and (4.24), the results in [13, Sects. 2, 4] show that also 
(-^0,n + A7q)~^ are integral operators whose integral kernels satisfy (4.25). □ 

Positivity improving and nondegeneracy of the ground state are considered next: 

Corollary 4.5. Under the hypotheses of Theorem 4.4, and under the assum,ptions 
that n is connected and Lei,J2 7^ Lq^^q, the positivity preserving relations in (4.21) 
a,nd (4.22) actually extend to positivity improving relations, thai is, =4 in (4.21) 
and (4.22) can be replaced by ^. In addition, the infimum of the spectrum of 
Ld,q., Lq.^q,, or Ln^u is a simple eigenvalue and the associated eigenfunction can 
be chosen to he strictly positive a.e. in fi. 

Proof. Positivity improving rather than just positivity preserving of the semigroups 
(resp., resolvents) in (4.21) (resp., (4.22)) follows from Corollary 2.14 and the fact 
that i^^Q is positivity improving if Ct is connected. The latter fact is implied, for 
instance, by the explicit lower bound for the Green's function GD,n{', •) = L^^-^{-, •) 
established in [48], [49], [140] for n ^ 3. Alternatively, one can invoke the 2>G 
theorem (proved, e.g., in [36, Theorems 6.5, 6.15], see also [3]. [37]) to obtain strict 
positivity of Gd,q,{', ■) onVlxQ,. In this context we also note; that indecomposability 
of e~*^°'" also follows from [50, Theorem 3.3.5] and positivity improving oie~^^'^-" 
and e~*^"'" is proved in [118] and [119, Theorem 4.5]. 

Nondegeneracy of the groundstate of Ld,q., Lq.^q, or L]v,o, and an associated 
strictly positive eigenfunction a.e. in fl then follows from [128, Theorem X. III. 44]. 
In this connect one recalls that the spectra of Lu.n, ^e^.n, and ^Af^n are purely 
discrete (cf. (3.59) and (3.67)). □ 

For background literature on nondegenerate groundstates we refer, for instance, 
to [46], [47, Ch. 7], [51, Ch. 13], [58, Sect. 10], [59], [67], [76, Sect. 3.3], [77], [128, 
Sect. XIII.12], [143, Sect. 10.5]. 

Corollary 4.6. Under the hypotheses of Theorem 4.4, all semigroups in (4.21) 
are sub-Markovian and hence extend to contraction semigroups on L°°{il;d"'x). 
Moreover, all semigroups in (4.21) extend to strongly continuous semigroups on 
LP{Q;d'^x), pe [l,oo). 

Proof. This is an immediate consequence of the interpolation and duality consid- 
erations discussed in [119, p. 56-57], upon noticing the following facts: 1) The 
generators of all semigroups in (4.21) are self-adjoint and nonnegative, and hence 
i^(0; c?"a;)-contractions. 2) Since e^*^™ " extends to an L°°(ri; d"a:)-contraction 
by [119, Corollary 4.10], the fact (2.83) then implies the L°°(n; d"a;)-contractivity 
of all remaining semigroups in (4.21). □ 
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Remark 4.7. One observes that condition (4.20) is automatically satisfied in the 
special case of local Robin boundary conditions considered in Lemma 3.11. In this 
context of local Robin boundary conditions a fair number of references proving 
Gaussian heat kernel bounds for L^i^n on bounded Lipschitz domains have been 
established in the literature as detailed in the paragraph preceding (1.2) in the 
introduction of this paper. The nonlocal Robin boundary conditions in terms of 
9 as encoded in (3.29) originated in [69] and to the best of our knowledge, the 
corresponding heat kernel and Green's function estimates for ^0,0 in Theorem 4.4 
are new. 

Remark 4.8. One can add a nonnegative potential ^ F G Ll^^{il; x) to all 
operators in Theorem 4.4 by employing the following standard procedure: First, 
adding the sesquilinear form £lv,n defined by 

£iv,a{u,v)= [ d''xV{x)u(^v{x), e dom (y^/^), (4.26) 

Jn 

to He.n and QD,n with domains H^{n) n dom (F^/^) and H^{fl) n dom{V^/^), 
respectively, yields densely defined and closed forms bounded from below. The 
uniquely associated nonnegative self-adjoint operators associated with 

£le,n{u,v)+Qv,Q{u,v), u,v e Hl{n) ndom{V^/'^), (4.27) 

and 

£1dMu,v)+QvMu,v), u,vGH\n)ndom{V^/^), (4.28) 

in obvious notation, will be denoted by Hqxi (and Hj^ q if = 0) and Hd q. 
One notes that C°°(j!) (the restrictions of C~(R")-functions to fl) and C^{n) 
are form cores for H@^q and HD,n, respectively. Temporarily replacing V by the 
bounded approximants = [1 — exp{—eV)]/e, e > yields positivity preserving 
semigroups e~*(^0'"+^=), e~*(-'^e,s2+^e)^ (,-t{L,^^a+Ve) ^ t ^ 0, employing the Trotter- 
Kato formula, 

g-t(A+B) ^ r^-t(A/m)^-t(B/m)im^ ^ ^ ^9) 

m— ^oo 

with ^ ^ = A* and B = B* E B(H). Using the monotone convergence for forms 
(cf. [28, Lemma 4] and the corresponding strong convergence of the semigroups 

g-t(iD.s!+i4)^ (^-t{LeM+v,) ^ g-t(Ljv,n+i4) e^*"'^-", e-*-H"e.n^ g-tH]v,n^ respectively, 

as e 4- (cf. e.g., [29, Lemma 5.2.13]), yields positivity preserving of the latter for 
all t ^ 0. (Alternatively, one can apply the truncation procedure for V described, 
e.g., in [58, Sects. 8, 9].) An application of Lemma 2.7 (Hi) then again yields the 
trace class property of all semigroups involved. Thus, the analogs of (4.21)-(4.22) 
hold for Hn,n, H@^^i, and i?jv,n- Moreover (cf. [46, Lemma 1.1], [28, Sect. 5.B]), the 
Trotter-Kato formula applied to form perturbations (either in the form established 
in [93] or using again approximations via Ve) also yields 

g-tHiv,n ^ g-tiw,fj^ t^O, (4.30) 

and hence by Theorem 2.3, the analogs of the integral kernel estimates (4.24) and 
(4.25) hold. 

One can also add a nonpositivc potential ^ W E £[^,^(^2; d^x) cither by apply- 
ing the perturbation method treated in [128, Theorem XIII.45] or by following the 
small form perturbation approach discussed in [58, Sect. 10]. This will, in general. 
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result in an additional multiplicative factor of the type e"^' on the right-hand side 
of (4.24). 

In this context of additive perturbations we also refer to [14] and [52], where 
Gaussian heat kernel bounds are obtained via Holder inequalities and domination 
of semigroups. 

Remark 4.9. For a variety of interesting implications of the i^-theory of Gaussian 
heat kernel bounds to L^-spectral theory and analyticity of semigroups we refer to 
[119, Ch. 7] and the extensive literature cited in the notes to this book chapter. 
Here we only mention the following facts (cf. also [23], [109], [110], [116], [119, p. 
95-97], [120], [136]): Introducing the sector S{4>) = {z e C\{0} | | arg(z)| ^ (j)}, (j) e 
(0,7r/2], all semigroups in (4.21) extend to holomorphic semigroups in LP{n;d^x), 
p e (1, oo), on the sector S{{n/2) - arctan [\p - 2|/[2(p - 1)^/^])). In addition, 

l|e"''""1lB(L.(n;.".)) ^ e 5((7r/2) - arctan (Ip- 2|/[2(p- 1)V2])), 

pe(l,oo), (4.31) 

and 

<y{Lp,n) C {0 e C I |arg(0)| ^ arctan (|p - 2|/[2(p - 1)^2]))} y {0}, p e (l,oo). 

(4.32) 

Here Lp^n denotes any of the generators of (the extension of) the semigroups in 
(4.21) on LP{Q.; (Tx), p e (1, 00). 

Finally, we mention a canonical counter example concerning positivity preserving 
semigroups among the set of nonnegative self-adjoint extensions of a strictly positive 

minimal operator: 

Remark 4.10. The literature on positivity preserving semigroups associated with 
the Friedrichs extension of differential operators is rather extensive as, typically, 
the Friedrichs extension corresponds to the case of Dirichlct boundary conditions. 
In the case where the associated minimal operator is nonnegative, it also possesses 
a second nonnegative distinguished self-adjoint extension, the Krein-von Neumann 
extension. The latter, in the context of Laplacians, was recently discussed in detail 
in [19] and [70] (see also [20] and the numerous references to the literature in 
these three sources) in the special case where the underlying minimal operator 
is in fact bounded from below by e > 0. Then, admittedly, on a formal level, 
the Krein-von Neumann extension in this special case appears to be defined in 
terms of a boundary condition that on the surface resembles a nonlocal Robin 
boundary condition, where O is expressed in terms of the operator-valued Dirichlet- 
to-Neumann map (cf. [19, Subsect. 5.2], [70, Sect. 13]). However, the well-known 
degeneracy of the lowest point in the spectrum, namely zero, of the Krein-von 
Neumann extension corresponding to the Laplacian associated with a bounded, 
connected domain f2 C M", n S N, and appropriate regularity of dU,, stemming 
from the nontrivial nullspace of the adjoint of the underlying minimal operator, 
guarantees that its semigroup is, in fact, never positivity preserving (let alone, 
positivity improving). For additional results in this direction, and the fact that 
Krein semigroups associated with L and bounded domains Q. are not Markovian, we 
refer to [66, Sect. 3.3]. The boundary conditions leading to a positivity preserving 
semigroup in the case of the (generalized) one-dimensional Laplacian on a bounded 
interval were classified in [60] (see also [66, p. 147]). 
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On an abstract level, and by appealing once more to [128, Theorem X. III. 44], 
the Krein von Neumann extension Ak of a densely defined and strictly positive 
operator A with deficiency indices d > 2 in some fixed complex, separable Hilbert 
space H will always have the rf-dimensional nuUspace of A* as its nullspace (cf., 
e.g., [19, Sect. 2], [20, Sect. 2]) and hence can never be positivity improving as long 
as the resolvent (or semigroup) of Ak is known to be irreducible. 



Appendix A. Sobolev Spaces on Lipschitz Domains in a Nutshell 

Following [69], we recall some basic facts in connection with Sobolev spaces on 
Lipschitz domains and on their boundaries. 

We start by briefly considering open subsets fl C M", n € N. For an arbitrary 
m e N U {0}, we follow the customary way of defining L^-Sobolev spaces of order 
±m in fl as 

H'^ifl) := {u e cTx) \ d^u e L^{Q.\ cTx) for ^ |a| < m}, (A.l) 

H-"'{n) := lu e V'{n) u= J2 '5"""' ^ith G L^iQid^'x), |a| ^ ml, 

(A.2) 

equipped with natural norms (cf., e.g., [2, Ch. 3], [111, Ch. 1]). Here V^fl) denotes 
the usual set of distributions on 17 C R". Then one sets 

H^(f7) := the closure of C^{n) in i?"(f7), m G N U {0}. (A.3) 

As is well-known, all three spaces above are Banach, reflexive and, in addition, 

(H^in))* = H-'^in). (A.4) 

Again, see, for instance, [2, Ch. 3], [111, Sect. 1.1.15]. Throughout this paper, we 
agree to use the adjoint (rather than the dual) space X* of a Banach space X. 

One recalls that an open, nonempty, bounded set fl C M" is called a bounded 
Lipschitz domain if the following property holds: There exists an open covering 
{Oj}i^j^N of the boundary dCl of Cl such that for every j G {1, . . . , TV}, Oj (1 fl 
coincides with the portion of Oj lying in the over-graph of a Lipschitz function (pj : 
M"~^ — ^ M (considered in a new system of coordinates obtained from the original 
one via a rigid motion). The number max{||V(/?j||ioo(j{,i-i.^n-i2./) 1 1 < j ^ N} is 
said to represent the Lipschitz character of f2. 

The classical theorem of Rademacher of almost everywhere differentiability of 
Lipschitz functions ensures that, for any Lipschitz domain O, the surface measure 
d^~^oj is well-defined on dil and that there exists an outward pointing normal 
vector p at almost every point of dQ. 

In the remainder of this appendix we shall assume that Hypothesis 3.1 holds, 
that is, we suppose that C K", n G N, n ^ 2, is a bounded Lipschitz domain. 

As regards i^-based Sobolev spaces of fractional order s G M, in a bounded 
Lipschitz domain c M" we set 

iJ^(M"):=|[/G5'(M-) ||f/||^,(Mn)=^^d"e|f/(0|'(l + ICI'1 <oo|, (A.5) 
H'{n) := {u G V'{n) \u = U\n for some U G if"(K")}. (A.6) 
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Here 5'(M") is the space of tempered distributions on M", and U denotes the 
Fourier transform of U G <S'(IR"'). These definitions are consistent with (A.l), 
(A. 2). Moreover, so is 

iJo(f^) {u € i?"(M") I supp(u) CH}, .s G R, (A.7) 

equipped with the natural norm induced by iJ'*(R"), in relation to (A. 3). One also 
has 

{H^{fl))* = H-%n), seR (A.8) 

(cf., e.g., [90]). For a bounded Lipschitz domain Q c M" it is known that 

{H'in))* = H-%fl), -l/2<s<l/2. (A.9) 

See [139] for this and other related properties. 

To discuss Sobolev spaces on the boundary of a Lipschitz domain, consider first 
the case when C M" is the domain lying above the graph of a Lipschitz function 
ip: M"-i R. In this setting, wc define the Sobolev space H^dn) for =^ s < 1, 
as the space of functions / G L'^{dfl: cf'^^Lu) with the property that f{x',(p(x')), 
as a function of x' E R"~^, belongs to i/*(M"^^). This definition is easily adapted 
to the case when il is a Lipschitz domain whose boundary is compact, by using a 
smooth partition of unity. Finally, for — 1 ^ s ^ 0, we set 

H'{dii) = {H-'idn))*, -l^s^O. (A.IO) 

Prom the above characterization of H" (d^) it follows that any property of Sobolev 
spaces (of order s £ [—1, 1]) defined in Euclidean domains, which are invariant under 
multiplication by smooth, compactly supported functions as well as compositions 
by bi-Lipschitz diff'eomorphisms, readily extends to the setting of H^{dn) (via 
localization and pull-back). As a concrete example, for each Lipschitz domain CI 
with compact boundary, one has 

H%dn) ^ L^idfl; dJ'-^u}) compactly if < s < 1. (A.ll) 

For additional background information in this context we refer, for instance, to [21], 
[22], [55, Chs. V, VI], [80, Ch. 1], [113, Ch. 3], [144, Sect. 1.4.2]. 

Assuming Hypothesis 3.1, we introduce the boundary trace operator 7^ (the 
Dirichlet trace) by 

^%:C{Ti)^C{dn), j%u = u\en. (A.12) 
Then there exists a bounded linear operator -yn 

7r,: H%n) H'-^^/'^\dn) {89,- (P-'^ u) , 1/2 < s < 3/2, 

7r, : H^''^{n) H^'^dQ) ^ L'^{dQ; dJ'-^uj), e e (0, 1) '^^''^^^ 

(cf., e.g., [113, Theorem 3.38]), whose action is compatible with that of 7^. That 
is, the two Dirichlet trace operators coincide on the intersection of their domains. 
Moreover, we recall that 

7r, : H'{n) H'-^^/'^\dn) is onto for 1/2 < s < 3/2. (A.14) 

Next, retain Hypothesis 3.1 and assume that 

AG^{H'{n)), l/2<s<3/2. (A.15) 

We then introduce the operator 7^^ (the strong Neumann trace) by 

^^ = u--fDAV: H'+\il)^L'^{dil;d''-'^uj), 1/2 < s < 3/2, (A.16) 
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where v denotes the outward pointing normal unit vector to d^. It follows from 
(A. 13) that 7jv is also a bounded operator. We seek to define the action of the 
Neumann trace operator in other (related) settings. Specifically, introduce the 
weak Neumann trace operator 

: {m e B^{^) I Lau e L^{^)) ^ H'^l'^idn), (A.17) 

as follows: Given u G H^{Q) with Lau £ L'^{Q) we set 

{<f>,lNu)y2 = + ($,LAw)L2(n), (A.18) 

for all (j) e ili/2(a0) and $ e H'^{n) such that jd^ = We note that the 
definition (A.18) is independent of the particular extension $ of (f), and that is 
bounded. 

Appendix B. Sesquilinear Forms and Associated Operators 

In this section we describe a few basic facts on sesquilinear forms and linear 
operators associated with them. Let "H be a complex separable Hilbert space with 
scalar product (• , • )-h (antilinear in the first and linear in the second argument), 
V a reflexive Banach space continuously and densely embedded into H. Then also 
H embeds continuously and densely into V*. That is, 

v^n^v*. (B.i) 

Here the continuous embedding T-L ^V* is accomplished via the identification 

■H3v^{-,v)n^V\ (B.2) 

and we use the convention in this manuscript that if X denotes a Banach space, 
X* denotes the adjoint space of continuous conjugate linear functionals on X, also 
known as the conjugate dual of X. 
In particular, if the sesquilinear form 

v(-, Vx V* -)-C (B.3) 

denotes the duality pairing between V and V*, then 

v{u, v)v. = {u, v)-H, u e V, V & U^V*, (B.4) 

that is, the V, V* pairing v( " ) " )v is compatible with the scalar product ( • , ■)n 
inn. 

Let T e B{V, V*). Since V is reflexive, (V*)* = V, one has 

T:Vh-V*, T*:V^V* (B.5) 

and 

v{u,Tv)v = v*{T*u,v)(v)' = v{T*u,v)v = v{v,T*u)v*- (B.6) 
Self-adjointness of T is then defined by T = T*, that is, 

v{% Tv)v' = V {Tu, v}v = v{v,Tu}v, u,vGV, (B.7) 

nonnegativity of T is defined by 

v{u,Tu)v ^0, ueV, (B.8) 

and boundedness from below ofT &y ct € M is defined by 

v{u,Tu)v, ^ ct\\u\\^, u€V. (B.9) 

(By (B.4), this is equivalent to v{u,Tu)v ^ ctv{u,u)v, u G V.) 
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Next, let the sesquilinear form a{ ■ , ■ ) : V x V C (antilinear in the first and 
linear in the second argument) be V-bounded, that is, there exists a Ca > such 
that 

\a{u,v)\^Ca\\u\\v\\v\\v, u,veV. (B.IO) 

Then A defined by 

A: ^ ~ ' (B.ll) 
\v>-^Av = a{-,v), 

satisfies 

A e B{V, V*) and v{u, Av)^, = a{u, v), u,vgV. (B.12) 
Assuming further that a( • , • ) is symmetric, that is, 



a{u,v) = a{v,u), u,v€:V, (B.13) 
and that a is V-coercive, that is, there exists a constant Co > such that 

a{u,u)^ Co\\u\\y, ugV, (B.14) 

respectively, then, 

A: V ^V* is bounded, self-adjoint, and boundedly invertible. (B.15) 
Moreover, denoting by A the part of ^4 in 'H defined by 

dom(A) = {u€V\Au€n} cn, A = A\^^^^^) : dom{A) H, (B.16) 
then A is a (possibly unbounded) self-adjoint operator in H satisfying 

A > Coin, (B.17) 
dom = V. (b.18) 

In particular, 

A-^ e B{H). (B.19) 

The facts (B.1)-(B.19) are a consequence of the Lax-Milgram theorem and the 
second representation theorem for symmetric sesquilinear forms. Details can be 
found, for instance, in [45, Sects. VI.3, VII.l], [55, Ch. IV], and [108]. 

Next, consider a symmetric form b{ ■ , ■ ) : V x V — )• C and assume that b is 
bounded from below by C}, € M, that is, 

6(u,m) ^ Cb||u||^, u€V. (B.20) 

Introducing the scalar product ( • , • )vi, : V x V — C (and the associated norm || • || y^) 

by 

{u,v)vi=b{u,v) + {l-Cb){u,v)H-, u,v£V, (B.21) 

turns V into a pre-Hilbert space (V; ( • , • )vb)j which we denote by Vfc. The form 
b is called closed in % if Vb is actually complete, and hence a Hilbert space. The 
form b is called closable in H if it has a closed extension. If b is closed in %, then 

\b{u,v) + {l-Cb){u,v)u\^\\u\\vM\v,, u,v€V, (B.22) 

and 

\b{u,u) + {l-Cb)\\u\\l\ = \\u\\l^, neV, (B.23) 
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show that the form b{ - , ■ ) + {1 — Cb){ - , ■ )'h is a symmetric, V-bounded, and V- 
coercive sesquilinear form. Hence, by (B.ll) and (B.12), there exists a Unear map 

[ vi-^ Bc^v = b{-,v) + {l-Cb){-,v)n, 

with 

Be, e B{Vb,Vt) and v,{u,B^,v)^, = b{u,v) + (1 - Cb){u,v)H, u,v&V. (B.25) 

Introducing the hnear map 

B = B,, + (C6 - 1)7: Vb ^ VI, (B.26) 

where I : Vb ^ VI denotes the continuous inclusion (embedding) map of Vb into 
V^, one obtains a self-adjoint operator i? in H by restricting B to 

dom(B) = {m € V I Bu e -H} C -H, B = B\^^^^s) '■ dom(B) H, (B.27) 
satisfying the following properties: 

B ^ Cbin, (B.28) 

dom (|B|i/2) ^ dom {{B - CblnY''') = V, (B.29) 

b{u,v) = {\B\^'^u,Ub\B\^'^v)^ (B.30) 

= {{B ^ CbluY'^u, [B - CbluY'^v)^ + Cb{u, v)-H (B.31) 

= v,{u,Bv)^,, u,v€V, (B.32) 

b(u, v) = (m, Bv)h, ueV, ve dom(B), (B.33) 
dom(i?) = {v & V\ there exists an fy G "H such that 

v) = (w, /„)« for all w e V}, (B.34) 

Bu = fu, u G dom(B), 

dom(B) is dense in H and in Vb- (B.35) 

Properties (B.34) and (B.35) uniquely determine B. Here Ub in (B.31) is the partial 
isometry in the polar decomposition of B, that is, 

B = Ub\B\, |B| = (B*B)i/2 > 0. (B.36) 

The operator B is called the operator associated with the form b. 

The facts (B.20)-(B.36) comprise the second representation theorem of sesquilin- 
ear forms (cf. [55, Sect. IV.2], [58, Sects. 1.2-1.5], and [94, Sect. VI.2.6]). 

Appendix C. On Heat Kernel and Green's Function Bounds 

In this appendix we briefly recall some bounds for heat kernels and Green's 
functions and prove the Green's function bounds as the latter in dimension n = 2 
are difhcult to find in the literature. 

In the case of Dirichlet boundary conditions (for general open sets f2 C M"), 

g-t(-AD,n) jg jjQown to be ultracontractive, that is, a bounded map from L'^i^l] d^x) 
into L°°{Q; d"a;), and to have a nonnegative integral kernel satisfying for all {x, y) e 
n X O, and ah f > 0, 

< e-*(-^°'"'(a;,y) < e-^""{x,y) = (47rt)-"/2e-l^-^l'/(^*' < {Ant)-"^^ (C.l) 
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due to domain monotonicity, as discussed in [50, Lemma 2.1.2, Example 2.1.8, 
Theorems 2.1.6 and 2.3.6]. Here Hq denotes the self-adjoint reaUzation of —A in 

i:2(M»;(i»a;), 

Ho = -A, dom(f/o) = H^{M."). (C.2) 
In addition, assuming Hypothesis 4.1, there exists a constant Cj^ao > Oi such that 
for all {x, y) G Cl X Q, and all 7 e (0, 1), t > 0, one has the bound 

< e-*^°'" {x, y) < c^,„„i-"/2e-l^-?'l V[4(i+7)ait] ^ (^.3) 

by [50, Corollary 3.2.8]. In particular, e~*^"'^ is ultracontractivc. 

Similarly, assuming again Hypothesis 4.1, also e"*^"-" is known to be ultra- 
contractive and to have a nonnegative integral kernel satisfying for some constant 
C7,ao,n > 0, all {x,y) e O x O, and all 7 e (0, 1), i > 0, 

< e-*^"'"(a;,2/) < C^,„o,n max(t-"/2, i)e-l^-2'lV[4(i+7)ait]^ (c.4) 

by [50, Theorems 1.3.9, 2.4.4, and 3.2.9] (see also the more recent [34], [15, Theorem 
4.4], [40], [118], [119, Ch. 6]). We note that the Lipschitz domain hypothesis is 
crucial in connection with the estimate (C.4) as it has to be modified in less regular 
domains characterized by Holder regularity of 90, as discussed in [31, Proposition 
3] (and the references cited in this context). 

For completeness, we also mention the explicit (but rather crude) Green's func- 
tion estimate in the context of the Dirichlet Laplacian and for Lo.n, based on 
domain monotonicity discussed, for instance, in [15], [49], [50, Ch. 3], [140] (see 
also [95, Sect. 1.2], [140]), and the references therein: For all {x,y) G fl x fl, and 
all A > 0, 

G%{-X,x,y) = {-AD,n + Mn)-\x,y) 



^ Go{-X,x,y) = {Ho + XIm.)-\x,y) 
2 

^(n-2)/2(A^/'k-t/|) 



r, I I X (2-n)/2 

2Tr\x — y\^ 
Ai/2 

^ jCx.n,n\x-y\^--, n^3, 

^ < X f^y, (C.5) 

\Cx,n\H^ + \x-y\-% n = 2, 

with Kv{-) the modified irregular Bessel function of order y (cf. [1, Sect. 9.6]). 
Similarly, in connection with LD,n one obtains, for all {x, y) gQ. x Q., and A > 0, 

CA,n|ln(l + |a;-y|-i)|, n = 2, 

The estimates (C.5) and (C.6) ignore all effects of the boundary dO. of O, but they 
suffice for the purpose at hand. 

Finally, we explicitly derive the analog of (C.6) in the case of Neumann boundary 
conditions as this is not so simple to find in the literature (particularly for n = 2). 

Lemma C.l. Assume Hypothesis 4.1. Then 

^ , , , fCao,oi,A,n,n|a;-t/|^"", n > 3, 

(-fAr,n(-A,a;,y) < 

lc„o,ai,A,n|ln(l + |ar-y|-i)|, n = 2, (C.7) 

A > 0, a;, y e O, x ^y. 
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More generally, let a G (0, {n/2)], then 

fao.a„aAO,„k-yP«-", aG(0,(n/2)), 
\Cao,aua,x,n\H^ + \x-y\-% « = n/2, (C.8) 

A > 0, X, y e fi, X ^ y. 

Proof. The main ingredients in deriving the bounds (C.7) and (C.8) are the identity 

(if + A/h)-" =r(a)-W dtt°'-^e-'^"+^\ a>0, A>0, (C.9) 

^0 

for a given self-adjoint operator fl" > 0, and the integral representation (see, e.g., 

[78, # 84327, p. 959]) 



f 

Jo 



dtt-^-^e-"*-'"^ ' = 2{h/a)-'/^K^{2{ahf'^), i/ € R, a > 0, 6 > 0. (C.IO) 

/o 

The asymptotic behavior of K^{-) (cf., [1, p. 375, 378]) 



- * ■Pii(i/2) + Ce][1 + 0{(i2)], i/ = 0, 



x4,o 



K,{x) = (7r/2)V2a;-i/2e-xri + o(/^-i-)i^ ^^0, (C.12) 
with Ce denoting Euler's constant (see, [1, p. 255]), then implies the bounds 

K„w4^['+^'^>^'"^""i['+'^'/"ri"'"'' .>o, (0.13) 

\c[ln(l+i--')][l + (2,T/i)>/2]-'e-, ^ = 0, 

for some universal constant C > in either case. Hence, combining (C.9) and 
(C.IO) with the heat kernel bound (C.4), and using the crude estimate, 

max(^-"/^ 1) < c^t-"/2e^*, e > 0, i > 0, (C.14) 

for some constant Ce > 0, then yields 

(7fjv,a + A7a)-"(x,y) 



pOO 

^0 

i^(n-2a)/2(((l + 7)«i)"'/'(A " e)'/^\x - y\) , 



^ a 



7,ao,ai,n,a,£: 



|a;_y|2l(2a-")/4 



A-£ 

' Cao,a^,a,X,n,e,n\x - t/P""", « € (0, (n/2)), 

^C'a„,ai,a,A,Q,£|ln(l + |a; - 2/1"^) |, Q! = n/2, 

A > £, a e (0, (n/2)], x,y€n, xj^y. (C.15) 

Here the last inequality follows from (C.13) and the fact that ft is bounded. Since 
£ > can be chosen arbitrarily small, this proves (C.7) and (C.8). □ 

We note again that the bound (C.7) for n ^ 3 can be found, for instance, in [54, 
Lemma 3.2] (see also [6], [35], [138]). 

We emphasize that the basic ingredients we used to derive (C.7) and (C.8) ap- 
pear, for instance, in [50, Sect. 3.4]. Moreover, the identical proof also yields the 
crude Dirichlet bound (C.6) (taking e = throughout). In our opinion, the value of 
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presenting Lemma C.l (besides establishing the bound for n = 2) hes in providing 
an explicit reference for the result and the ease with which it is derived. 
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